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1 Introduction 

It was Planck who first demonstrated the power of statistical methods in quantitatively under- 
standing a gauge theory pQ. His important suggestion was to subject indeterministic phase and 
amplitude changes of a single resonator in the wall of a cavity - in thermal equilibrium with the 
contained electromagnetic radiation - to an averaging procedure dictated by the laws of (sta- 
tistical) thermodynamics. Appealing to a known, classically derived result on gross features of 
the so-called black-body spectrum (Wien's displacement law), postulating a partitioning of the 
total energy into multiples of a smallest unit, and appealing to Boltzmann's statistical definition 
of entropy, Planck deduced his famous radiation law. As an aside, he discovered a universal 
quantum of action needed to relate the entropy (disorder) and mean energy of a single resonator 
to its frequency. The robustness of his result is demonstrated by the fact that even for physical 
objects sizably deviating from ideal black bodies Planck's radiation law holds to a high degree 
of accuracy. 

The purpose of the present article is to give a concise presentation of results, accumulated 
over the last four years, on generalizations of the thermal U(l) gauge theory studied by Planck: 
SU(2) and SU(3) Yang-Mills thermodynamics. It is possible that an SU(2) gauge symmetry, 
dynamically broken down to U(l) by a deconfining thermal ground state, underlies photon 
propagation [21 [3j HI [6] . We do not here consider SU(N) gauge theories with oo > N > 4 the 
reasoro being nonunique phase diagrams [3 [8]. 

Yang-Mills thermodynamics strongly relates to the concept of emergent phenomena. On the 
most basic level, temperature itself is an emergent phenomenon depending on the fluctuating 

x The case N= oo may be an important exception [7|- The nonuniqueness of the phase diagram is due to an 
incomplete Abelianization of the fundamental gauge symmetry by an adjoint Higgs field: There is no principle 
which decides at what temperature a nonabelian factor is broken to its Abelian subgroup. 
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degrees of freedom defining it. Conversely, the mass of a magnetic monopole, which, as a short- 
lived field configuration contributes to the thermodynamics of the Yang-Mills ground state at 
high temperatures, is determined by temperature. That is, a single monopole owes its existence 
to the existence of all other fluctuating monopoles and antimonopoles in the ensemble. 

In Yang-Mills thermodynamics various temperature-dependent emergent phenomena, facili- 
tated by topologically nontrivial mappings from submanifolds of four dimensional spacetime into 
the (sub)manifold(s) of the gauge group, dominate the ground-state thermodynamics in three 
different phases. Albeit their microscopic dynamics is complex and not accessible to analytic 
treatment, a thermodynamically implied spatial coarse-graining down to a uniquely determined 
resolution, determined by the Yang-Mills scale and temperature, yields accurate and technically 
managable representations of the (ultraviolet-regulated) partition function at any given tempera- 
ture@. Here the term spatial coarse-graining refers to the process of eliminating short- wavelength 
gauge-field fluctuations in favor for effective fields and their couplings in a reformulation of the 
same partition function at lower and lower resolution and at a given temperature. This coarse- 
graining leads to the emergence of an effective action S eS being a functional of the effective fields 
which determines the weight for the functional integration over the latter in the reformulated 
partition function valid at a given maximal resolution. 

Remarkably, starting out from exact BPS saturated solutions to the euclidean field equations 
in the deconfining phase, the derivation of the thermal ground state in the deconfining, high- 
temperature phase makes no reference to the way of how the continuum partition function of 
the theory is regularized in the ultraviolet. This is a consequence of the possibility to compute 
the expectation value of a uniquely determined operator, representing the phase (ft of an adjoint 
scalar field (ft, over suitable, topologically nontrivial, and BPS saturated configurations first 
and to subsequently investigate its average effect on integrated fluctuations and its direct effect 
on explicit field configurations in the topologically trivial sector. No reference is made to an 
ultraviolet regularization in the process. When decreasing the maximal resolution in the BPS 
saturated situation it is possible to deduce definite information about eft's modulus \cj>\ from its 
phase (ft since the process of spatial coarse-graining, which for (ft rapidly saturates the limit of 
vanishing maximal resolution at a finite maximal resolution, yields a spacetime homogeneous 
(and of course isotropic) valu^l of \(ft\ starting from a maximal resolution determined by \(ft\. 

That is, spatial coarse-graining over noninteracting, BPS saturated configurations in the (only 
admissible) sector with topological charge modulus \Q\ = 1 dimensionally formally reduces this 
sector of the theory from D = 4 to D = 1 (quantum mechanics with periodic-in-time trajecto- 
ries) at a finite maximal resolution. Thermodynamically, the gauge- invariant quantity \(ft\ must 
not carry energy (or momentum in 4D) because of its inherited BPS saturation. This is also the 
reason why no local vertices of the field (ft with the topologically trivial, coarse-grained gauge 
field involving three of more external legs of the latter may exist: Such vertices would on the 
level of the effective theory convey energy-momentum transfer from the topologically trivial to 
the topologically nontrivial but BPS saturated sector of the theory which contradicts the very 
existence of the field (ft. Notice, however, that on the fundamental level topologically trivial 
fluctuations do interact with the topologically nontrivial sector exchanging energy-momentum 
associated with a resolution higher than \(ft\ thus having no visible effect on the field (ft. These 
interactions, introducing a temporary (anti)caloron holonomy in addition to short-range radia- 
tive corrections, are described by a pure-gauge configuration in the effective theory which in fact 
lifts the energy density of the thermal ground-state from zero to a finite positive value. 

Despite the fact that the Stefan-Boltzmann limit is approached in a power-like and thus rapid 

2 An exception is the nonthermal behavior shortly below the Hagedorn transition. 

3 Due to the nonfluctuating nature of the field (j> or higher Lorentz-spin fields, potentially generated by a coarse- 
graining over the BPS sector of the fundamental theory, infinite-volume thermodynamics excludes the existence 
of the latter and demands homogeneity (constancy) of the former's gauge invariant modulus. 
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way with increasing temperature the inherent resolution \<p\ of the thermal Yang-Mills system 
decreases with temperature, and thus any memory of the short-distance regularization of the 
partition function is wiped out, see also [9] for the corresponding lattice observation (increasing 
derealization of topological charge with increasing temperature). Conceptually, this is in accord 
with the situation known in zero-temperature perturbation theory where our ignorance about the 
values and UV regularization dependence of bare parameters is shown to be no obstacle to the 
predictivity of the quantum Yang-Mills theory since, upon their dressing at a finite resolution, 
only finitely many such parameters need to be fixed |1Q|. [TTT [12], H3] . 

Following the appreciated advice of a Referee, the presentation in this paper resorts to a 
hybrid style: A statement, whose validity is argued for in a more physical rather than rigorous 
mathematical way, and a number of definitions are highlighted by slanted script and are inter- 
spersed into the argumentation, a statement whose validity under stated assumptions is verified 
by direct calculation is presented according to custom within the mathematics literature, apolo- 
gies to the irritated reader. But even though the subject presented (4D quantum field theory) 
still awaits its rigorous mathematical foundation the pragmatic line of pursuit followed in the 
present work unearthens a number of unexpected, quantitatively very precisely representable 
facts which, as the author is convinced of, do not depend on a future, rigorous formulation of 
quantum field theory. 

For maximal benefit it is recommended to read the present article in conjunction with [7\. 
That article contains extended discussions of the involved physics, explicit expressions for the 
relevant topological field configurations, and graphic displays of numerical results. Although 
the use of differential forms would simplify certain statements in the beginning of Sec. [3] our 
presentation entirely resorts to the component notation. 

The following conventions will be used: Einstein summation (summation over doubly oc- 
curring indices), solely lower-case indices for contractions in the euclidean formulation, lower- 
case and upper-case indices for contractions in the real-time formulation, and natural units 
(h = ks = c = 1). 

The outline of this work is as follows: In Sec. [2] we remind the reader of basic facts about ther- 
mal Yang-Mills gauge field theory. Sec. [3] discusses the deconfining phase where a thermal ground 
state emerges upon a spatial coarse-graining over interacting, BPS saturated field configurations 
and where the Yang-Mills scale occurs as a purely nonperturbative constant of integration. We 
give tight estimates on the goodness of the finite-volume saturation of the coarse-graining pro- 
cess, and we account for the radiative corrections in the effective theory. The thermodynamics 
of the intermediate, preconfining phase is addressed in Sec.[H Here the unbroken abelian gauge 
symmetry of the deconfining phase is dynamically broken by monopole-antimonopole conden- 
sate^). Emphasis is put on a discussion of supercooling which takes place because the switch 
from small to large caloron/anticaloron holonomy is energetically disfavored. Finally, in Sec. [5] we 
elucidate the process of the decay of the monopole-antimonopole condensate of the preconfining 
phase giving rise to a zero-pressure and zero-energy density ground state in the confining phase. 
In that phase no propagating gauge modes exist, and the spectrum is represented by single or 
selfintersecting center- vortex loops which we interprete as spin-1/2 fermions. The naive series 
for the thermodynamic pressure represents an asymptotic expansion, and we show its Borel 
summability for complex values of the expansion parameter. Upon continuation to the physical 
regime a sign-indefinite imaginary part is encountered which for sufficiently small temperatures, 
however, is, smaller than the definite real part. 
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2 Thermal Yang-Mills theory 

2.1 Euclidean formulation and symmetries 

On a fiat, four-dimensional euclidean spacetime with coordinates < r < j3 = T _1 (time) and 
x (infinite three-dimensional space) the partition function Z of a pure Yang-Mills gauge-held 
theory subject to the gauge group SU(N) is formally dehned as 

Z= I VA^eM-S], (1) 

J j 4 m (t=0,x)= j 4 m (t=/3,x) 

where the gauge-held conhguration A^ is Lie-algebra valued, A^ = t a , (a = 1, • • • , JV 2 — 1), 
with the generators t a in the fundamental representation normalized as trt a tb = ^S a b , and T 
is the temperature. The action S is dehned as S = tr Jq dr J d 3 x F^ U F^ U where g is a 
dimensionless coupling constant, F^ v = d^A v — d v A^ — i[A^, A v \ , (\x,v = 1, • • • ,4), and the 
measure of (functional) integration is VA^ = X\ T g a dA a ^(T,x). Here the product is over the 
continuously varying values of the coordinates t, x, and over a = 1, • • • , N 2 — 1 

For the gauge group SU(2) we set t a = \\ a where A a are the Pauli matrices. The integration 
measure VA^ = \\ T ^dA^(T,x) is ill-defined as it stands. If the field theory is endowed with 
an ultraviolet regularization then the infinite product is over a discrete index. We will argue 
that Yang- Mills thermodynamics, formally defined @ by Eq. ([1]), does not relate to the way how 
sense is made of the formal object in Eq. ([I]) by introducing a minimal length as long as the 
selfadjusting resolution in the system is lower than the resolution associated with this ultraviolet 
scale. As already discussed in the Introduction at a given temperature T a unique maximal 
resolution \(j)\ appears to emerge. Since \<p\ decays like a power when increasing T it is guaranteed 
that the high-temperature physics is insensitive to any definite ultraviolet substantialization of 
Eq. (H]). The action density ^ tr F^ V F^ V in Eq. ([1]) is invariant under gauge transformations 

A^ ——* VlA^rfi + iQd^fl^ , where £1 is an element of SU(N) in the fundamental representation, 
but the functional integration is carried out over gauge- inequivalent, periodic- in-r gauge-field 
configurations. 

In lattice definitions of the partition function in Eq. ([T} one can show its invariance under 
temporally local center transformations, that is, under gauge transformations which are periodic 
up to a multiplication with a center element: 0(r = 0,x) = Z 0,(t = f3,x) where Z E Z N and 
thus the relevant gauge group actually is SU(N)/Z N . If this (electric) center symmetry of the 
Yang- Mills action S, subjected to the gauge group SU(N), is broken dynamically then a Z N 
degeneracy of the ground state must be present. Furthermore, the action S is invariant under 
continuous spatial rotations and translations. There is also an invariance w.r.t. time translations 
r -» r + r , < t < (3. 

2.2 BPS saturated field configurations at finite temperature 

The Euler-Lagrange equations, D^F^ V = (stationarity of the action, j^- = 0) are solved 

by configurations A^ obeying the (anti)selfduality condition F^ u = rfci^,. Here D^- = ■ 
— i[An, ■] and F^ = ^e^^F^x where e^ VK \ is the totally antisymmetric tensor with €1234 = 1. 
(Anti)selfdual configurations saturate the Bogomol'nyi bound on the action (BPS saturation): 
S = ^pr\Q\ where Q = Jq dr J d 3 x F® U F® U G Z is the topological charge. Q is finite and 
quantized if according boundary conditions are imposed. At finite temperature we consider BPS 
saturated, finite-action configurations A^ which behave accordingly at spatial infinity and are 



4 This partition function implies Legendre transformations between formally defined thermodynamical quanti- 
ties like pressure (minus free energy), energy density, and entropy density. 
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periodic in r. In this case Q is an integer. As is usual, we refer to these configurations as 
calorons/anticalorons for a positive/negative sign of Q. (The presence of spatial boundaries, 
which are not at infinity, would explicitely break the translational invariance. We so far have 
no well backed up insight on how to treat Yang-Mills thermodynamics analytically in this case 
which, in a softend fashion, is also conveyed by external sources. A possibility for the treatment 
of mild distortions is to deform the undistorted situation adiabatically by letting T — > T(x) 
and/or A — ► A (if) where A denotes the Yang-Mills scale). Solutions to the Euler-Lagrange 
equations, which are not (anti)selfdual, have been constructed numerically, for recent work on 
axially symmetric finite-temperature and instanton/antiinstanton configurations and see |14j 
and [15], respectively. 

Proposition 1. The euclidean energy-momentum tensor 9^ v = —F^ X F^ X + F%\F%^ and 
every local, scalar composite of the form trt a F^F^, trt a F^ K F KL/ F U ^, trt a F^F^F^Fx^, ■ ■ ■ 
vanish identically on a BPS saturated field configuration A^. 

Proof. Routine computation. ■ 
Notice that for Q = we have S = implying that calorons in this sector are pure gauges: 

The Polyakov loop P(x)[A] is defined as P(x)[A] = V exp[i A±(t, x)] where the symbol V 
demands path-ordering. 

A caloron/anticaloron is said to be of trivial holonomy if 

Poo[A] = limi^i^oo P(af)L4] G center of the gauge group. 

(Recall that for SU(N): center={exp[2m^]l N \k = 0, 1, • • • , N- 1} = Z N .) 

Notice that for the gauge group SU(N) this definition does not depend on the choice of gauge as 
long as n(r = 0,x) = ZQ(t = (3,x), where Z £ Z N , since P^A] Q(t = 0)P oo [A]^(t = (3) . 

Example 1. For the gauge group SU(2) the following calorons/anticalorons (Harrington-Shepard 
(HS) [16]) are of trivial holonomy and of topological charge Q = ±1: 

A^(t, x) = rf t a d y lnll(r, r) (caloron, Q = +1) , 



A^(t, x) = rf^ v t a d v lnn(r, r) (anticaloron, Q = —1) , (2) 

where the 't Hooft symbols n afMU and 77* are defined as r/" = e a ^ u + 5^5^ — 5% 6^4 and ff^ v = 
e fiu ~ ^VvA + <^"^4; an d the prepotential II is given as 



2 shih C^nr") 

H(r,r)^l + ^ X / 9 N » ( r = 1^1) • 

Pr cosh hnr) - COS (?f) 



The dimensionful modulus p is inherited from the singular-gauge instanton configuration with 

2 

prepotential lip (r, r) = 1 + T i + gi since IT is obtained from IIo by superimposing its infinitely 
many images in mirrors placed at r = and r = (3 to generate periodicity in r. Additional 
moduli are the shifts r — ► r + t z , (0 < t z < (3) , and x — > x + z and, if one wishes, global gauge 
transformations. 

Example 2. For the gauge group SU(2) there exist [T71 [181 (13 US [23] explicitely constructed 
|24 [ 1251 [26] \27 \ 128] BPS calorons/anticalorons (Lee-Lu-Kraan-van-Baal (LLKvB)) of nontrivial 
holonomy and topological charge Q = ±1. For the trivial- holonomy case configurations with 
\Q\ > 1 were constructed in [2T| [22] . 
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In contrast to their trivial-holonomy counterparts |29j these configurations are in isolation 
instable under quantum deformation [42J: For a holonomy sufficiently close to trivial and for 
p > the static BPS magnetic monopole and antimonopole constituents [39 [ I40 [ [4"T] attract under 
the influence of quantum fluctuations and thus eventually annihilate one another. This relaxes 
the LLKvB caloron or anticaloron back to the stable situation of a HS caloron or anticaloron. For 
a holonomy far from trivial and for p > BPS magnetic monopole and antimonopole repulse [42] 
one another under the influence of quantum fluctuations. As a consequence, the large-holonomy 
LLKvB caloron or anticaloron dissociates into a pair of a screened BPS magnetic monopole and 
its antimonopole. Screening occurs due to the presence of short-lived magnetic dipoles that are 
provided by intermediary small- holonomy LLKvB calorons and anticalorons. 

In the seminal work [29] nontrivial- holonomy calorons/ anticalorons were argued to not con- 
tribute to the partition function based on the observation that quantum corrections produce a 
term in their effective action which diverges like the three-volume of the system. This argument 
is certainly correct if holonomy is considered a quantity that is externally sustained at a definite 
value. Hence no explicit integration over the holonomy must occur in any first-principle evalua- 
tion of the (ultraviolet regularized) partition function at sufficiently large temperature. However, 
viewed as a dynamical, short-lived quantity nontrivial holonomy does occur through the quan- 
tum induced deformation of the trivial-holonomy case. Partially based on theoretical work on 
fermionic zero modes, used as a diagnostics for lumps of topological charge [HUSH] to avoid the 
application of a cooling procedure to a given configuration, this is impressively demonstrated 
by many lattice investigations. A nonexhaustive list of references is [U [30l [HJ E2J I33] [34] . 

The construction of nontrivial- holonomy calorons/anticalorons of higher topological charge 
for a Yang-Mills theory subject to the gauge group SU(N) was investigated by Bruckmann and 
van Baal in [37J and the explicit form of the solution was given for \Q\ = 2. The interesting result 
is that these calorons possess n\Q\ constituents monopoles whose sum of magnetic charges (with 
respect to U(1) N " 1 ) is nil. Recently, a nontrivial-holonomy caloron of Q = 2 and nonvanishing, 
overall magnetic charge was constructed [38]. 

Remark 1. On the level of the euclidean saddlepoint one has for the masses mi of a BPS 
magnetic monopole and m,2 of its antimonopole inside a LLKvB caloron: m\ + rri2 = 8ir 2 T 
[24] . Thus already on the classical level one observes the remarkable fact that the emergence 
of a particular monopole depends on the emergence of temperature or in other words on the 
existence of all other fluctuating monopoles and antimonopoles and propagating gauge fields in 
the ensemble. 



2.3 Propagating fields at finite temperature: Q=0 

In a given gauge and in the euclidean formulation the topologically trivial sector {<L4^} is 
represented by a superposition of plane waves: 



5A^(t,x) = ^2 exp 



2irin 



P 



6Au, n (x) 



(3) 



where 8A^ n {x) = f d 3 ka^ n (k) exp[i/c • x] and the function a Mi „(/c) falls off sufficiently fast in 
\k\ and in n for the integrals and the sum in Eq. (|3|) to exist, respectively. The quantity is 
called nth Matsubara frequency. 



Remark 2. Upon a Wick rotation r 



it . 



(t real) , one shows that the propagator of the 



field SA^ decomposes into a quantum part (describing a particle of four-momentum p possibly 
being off the mass shell, p 2 = p^p^ ^ 0) and a thermal part (describing thermalized on-shell 
propagation), see for example [33]. 
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3 Deconfining phase 

3.1 Thermal ground state: Interacting calorons/anticalorons 

If not stated otherwise we consider the gauge group SU(2) from now on. We perform a spatial 
coarse- graining over the sector of nontrivially BPS saturated (nonpropagating) field configura- 
tions in singular gauge to arrive at a nonpropagating adjoint scalar field 4> of spacetime indepen- 
dent modulus. Our strategy is to derive (f>'s second-order equation of motion and, by requiring 
compatibility with BPS saturation, to subsequently determine the field (f> (modulus and phase) 
in terms of T and a constant of integration A. The perturbative renormalizabilit}0 of the sector 
with propagating gauge fields (Q = 0) [lOl [111 021 H3] and the requirement of gauge invariance 
then yield a unique effective action which is associated with a maximal resolution given by 0's 
modulus. 

If in the effective action a spatially homogeneous composite held emerges after spatial coarse- 
graining over the sector of nontrivially BPS saturated field conhgurations of trivial holonomy 
then this composite is a scalar under rotations (0(3) scalar) and transforms in the adjoint 
representation of the gauge group SU(2). 

Note 1. The ^-component and (nonlocal) products thereof are 0(3) scalars only in covariant 
gauges. 

As stated in the Introduction, the term 'spatial coarse-graining' refers to a lowering of the 
maximal resolution available in the system at a given temperature when keeping the partition 
function fixed. The term 'effective action' refers to minus the exponent in the weight accord- 
ing to which an average over configurations is performed in the partition function after spatial 
coarse- graining. 

exclusion of explicit nontrivial holonomy: As shown in |29j . calorons/anticalorons with explicit 
nontrivial holonomy induce a one-loop effective action which diverges with the three-volume of 
the system. Thus explicit holonomy must not enter the process of spatial coarse-graining in the 
caloron/anticaloron sector. (Short-lived implicit holonomy, however, emerges by quantum defor- 
mation of the trivial-holonomt case and is responsible for the generation of short-lived magnetic 
dipoles (large thermodynamic weight) or screened magnetic monopoles and antimonopoles (very 
small thermodynamic weight). 

transformation property under 0(3): Since nontrivially BPS saturated field configurations of 
trivial holonomy are nonpropagating field configurations their coarse-grained counterparts rep- 
resent spatially homogeneous background fields in the effective theory. But the existence of 
a nontrivial 0(3) tensor after coarse-graining would spontaneously break rotational invariance 
which is impossible in the absense of microscopic degrees of freedom that single out a direction 
in space at vanishing momentum. 

gauge transformation property: The scalar <f) must transform homogeneously under a change of 
gauge for otherwise the coarse-grained gauge field 5A^ would have to form a composite to couple 
to <p in a gauge-invariant way. The existence of such a composite on the level of the effective 
action would, however, contradict perturbative renormalizability |10| [TT1 [T2| [13] which states 
that all propagating degrees of freedom are represented by 5A^ itself. The only homogeneously 

5 This is the statement that the only effect of integrating out fluctuations in the sector with Q = is a resolution 
dependence of the gauge coupling and the normalization of a plane wave in a given gauge. Modulo these radiatively 
generated effects the effective action describing fluctuations in the sector with Q = has the same form as the 
fundamental action. In the absence of external sources probing the thermal system its inherent resolution \<f>\ is 
a function of temperature. As a consequence, the effective gauge coupling is a function of temperature and, in 
the (only) physical gauge, the wave function normalization is the characteristic function \: \ equals unity if the 
plane wave resolves its environment by less than \4>\, and \ equals zero if the plane wave would resolve by more 
than \(f>\ since such a fluctuation already is integrated out. 
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transforming, nontrivial quantities in the fundamental theory are (nonlocal) products of the field 
strength F^. Since 

tatb = ^{ta,t b } + ^[t a ,t b ] = ^ Q<5 afe l 2 + ^afec^ (4) 

we may without restriction of generality schematically write 

f 1 -* ~ tr (t ai ■ ■ ■ t aK F ■ ■ ■ F ) , (K > 1) 

with appropriate contractions of the Lorentz indices and parallel transports for the field strength 
F = F b tb implied. By virtue of Eq. ([4]) this can always be decomposed into spin-0 and spin-1 
representations of SU(2). The case of the spin-0 representation (gauge- invariant composite) is 
irrelevant because it decouples and no energy- momentum is associated with it, see proposition 
HJ Thus we are left with the spin-1 representation which proves the claim. 



Notation: We denote by {(r, 0), (r, x)} the spacelike Wilson line V exp i J^q\ dz^A^^z) 
the path of integration is a straight line. 



where 



The following is the unique definition for the set K, of r-dependent algebra-valued functions 
which contains (ft's phase (ft = A- , (|0| 2 = tr \ (ft 2 ): 

K = { Yl j d * X j dp tr tF^fafyc {(r,0),(r,f)} Q F^(r,f) Q {(r, x), (r, 0)}^ } , (5) 

where the sum is over a HS caloron and anticaloron (singular gauge) and t = (i^t^i 3 ). 
Explicitly, the set /C is parametrized by the coordinates of the caloron/ anticaloron center 
zc,A = ( T C,A, zc,a)- The integrals are over infinite space and the entire range < p < oo 
of the scale parameter p, and, as we shall see, the set fC is implicitly parametrized by arbitrary 
rescalings and global gauge transformations. 

HS caloron/ anticaloron: 

Again, explicit holonomy is excluded by the result of the semiclassical calculation in [29] but im- 
plicit, short-lived holonomy emerges by the quantum deformation of the trivial-holonomy case. 
local dehnition: 

This is excluded by Proposition Q3 

curved path for evaluation of Wilson line {(r, 0), (r, x)}: 

Since the path is purely spacelike there exists no mass scale on the level of BPS saturated field 
configurations to parameterize curvature. 
higher n-point functions: 

Since JC contains the dimensionless phase (ft all its members must be dimensionless. Considering 
nonlocal, n-fold products of Fu V with n > 2 together with the associated additional space inte- 
grations, a factor of 1 ~ n would have to be introduced to make these contributions dimensionless. 
Since there are no explicit dependences on (3 on the level of BPS saturated field configurations, 
see beginning of Sec. 2.2, this possibility does not exist. 
moduli-space average: 

(i) Integrations over the dimensionful moduli p and tc,a, zc,a must have a flat measure since the 
members of KL make no reference to any scale on the level of BPS saturated field configurations. 

(ii) The right-hand side of Eq. ([5]) transforms in the adjoint representation. Shifting the calo- 
ron/anticaloron spatial center from to zq,a an d honoring spherical symmetry, an additional 
pair of Wilson lines would have to be introduced to parallel transport F^ u from (0, 0) to (0, zc,a)- 
Integrating then over zc A yields zero. (If this integral would not vanish then the scale (3 would 
occur explicitely in the definition of the dimensionless members of KL. But this is forbidden on 
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the level of BPS saturation, see beginning of Sec. 2.2.) 

(hi) An integration over tc,a yields zero. (The case of a constant term in the Fourier series asso- 
ciated with the integrand again would imply that the scale (3 occurs explicitely in the definition 
of the dimensionless members of /C.) 

(iv) Since the members of K, are gauge-variant objects integrations over the global gauge orien- 
tations of the caloron/anticaloron yield zero and thus are forbidden. 

shift — > y ^ 0: Shifting — > y in Eq. (JS} but leaving the caloron/anticaloron center fixed at 
0, spherical symmetry would imply the need for an additional pair of Wilson lines to connect 
y with 0. But this is just a global gauge rotation of the unshifted situation and thus does not 
alter /C. 

caloron/anticaloron with \Q\ > 1: 

Besides the translational moduli there are m > 1 dimensionful moduli in such a configuration. 
For example, a trivial-holonomy caloron with Q = 2 has three dimensionful moduli: two scale 
parameters and the distance between the two centers of its topological charge. Considering an 
n-point function (n nonlocal factors of the field strength F^ v ) with n — 1 integrations over space 
and a flat-measure integration over the m dimensionful moduli (not counting the shift moduli) 
of the caloron, we arrive at a mass dimension 2n — 3(n — 1) — m = 3 — n — moi the object. To 
avoid the introduction of explicit powers of (5 (BPS saturation) in the definition of JC this mass 
dimension needs to vanish. But for n > 2 and m > 1 we have 3 — n — m ^ 0. 



Proposition 2. The Wilson line |(r, 0), (r, x)\ 



C,A 



evaluates to 



|(r,?),(r,ic)| = cosg±2it b ^- sing, (6) 

where g = g{r,r,p) = g(/3f,/3f,/3p) = g{f,f,p) = J* ds §<9 T logII(r, sr, p). The + or - sign 
relates to a caloron or an anticaloron, respectively. Explicitly, one has 



g = —ir 2 p 2 sin(27rf ) / ds 
Jo 



1 1 sinh(2-7rfs) 



s [cosh(27rfs) — cos(27rf )] [cosh(27rfs) — cos(27rf ) + sinh(27rfs)] 

(7) 

The function g exists and approaches constancy in f more than exponentially fast with increasing 
f > 1. 

Note 2. To point out essential properties only we have set the phases tq and ta of the r 
dependences of caloron and anticaloron equal to zero. They can easily be reinstated by letting 

T -> T + TC,A- 

Proof. 

irrelevance of path-ordering: 

Observe that J^'^ dz fM A fM (z)\ c A = ± ds XiAiir, sx) = ±tbX b d T Jq ds log n(r, sr, p). That 

is, the integrand in the exponent of { (r, 0), (r, x) > varies along a fixed direction in the Lie 

I J c,A 

algebra of SU(2). Path-ordering thus can be omitted. 

explicit form of < (r, 0), (r, x) > : Routine computation. 

I- J C,A 

existence of g for all values of its arguments: 

The potentially problematic point in the domain of integration is s = for f = k E 7L. By 
Taylor expanding the sine function in front of the integral in Eq. ([7]) about f = k and by Taylor 
expanding the cosine and the cosine hyperbolic functions in the denominator of the integrand 
about f = k and s = 0, respectively, one easily checks that the limit f — * k exists when p > 
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and f > 0. 

saturation property for growing r > : 

Split the integration in Eq. © as 2 = f} ds = X\ + J2 = /o 2 *^ ^ s + ds = L dz + J x m dz. 
Z\ does not depend on f. The integrand I is given as 

z[cosh.z — cos(2-7rf )][cosh z — cos(27rf) H — ^p- sinhz] 
For the integration in I2 the integrand / is bounded from above as 

2 e z 

71 — 5 , 9 . Vf,p;z>l. (9) 

(e 2 — 2 cos(27rrjj z 

Since f^ r dz I(z, p, f) = dz I(z, p, f) — f^° f dz I(z, p, f) and since, by virtue of Eq. ([9]), the 
modulus of the second summand is bounded by *=* 2 — ttt-^t we are assured a more than 

e ^irr_2 COs(27rr) 

exponentially fast saturation in f. Numerically, ^ 2lTi ' ^fl^iMl < 10~ 5 for f > 2 independently of 
p. (Saturation in p is now trivial.) ■ 

Proposition 3. The integrand in Eq. |5j) ; when evaluated on a caloron, is as 



_ 2 32vr 4 x a vr 2 /5 4 + P 2 (2 + cos(27rf )) 
3 r (27r 2 p 2 + 1 - cos(2vrf)) 2 



■i^T /J^^ 2 >< F im , (io) 



where the functional F is given as 

Flm _ 2c « (25) fM - +8in(2ffl ( 2 JW _ + 



n 2 n y v y V n 2 n 2 n n 

(ii) 

Proof. Lengthy routine computation, see |44^ I45j . ■ 

Proposition 4. The integrand in Eq. ([5|) ; when evaluated on an anticaloron, is obtained by a 
parity transformation, x — > —x, of the integrand evaluated on a caloron. 

Note 3. For equal temporal phases, tc = ta, it then follows that the integrands cancel. But, 
as we will show, nontrivial BPS saturation of the field <ft requires that tq — ta = db5. 

Proof. It is easily checked that F ixv {t,x)c = F^ u (t, —x)a and that 

{(r,0),(r,f)} c = ({(r,f),(r,0)}J t = {(r, 0), (r, -x)} A = ({(r, -x), (r, 0)} J*. This 
proves the claim. ■ 

Remark 3. Due to the appearance of the factor — in the expression (|10p the unconstrained 
angular integration in Eq. (JH) yields zero. Thus for the final integration over p to possess a 
nonvanishing integrand the radial integral must diverge. 

Proposition 5. The only term in F[<7,II], which gives rise to the divergence of the radial 
integral, is — sin(2^)^J^. This divergence is logarithmic. 

Note 4. Since only spatial derivatives are involved this term arises from magnetic-magnetic 
correlations. But it is the magnetic sector whose insufficient screening gives rise to the poor 
convergence properties in thermal perturbation theory |47j . 
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Proof. Obviously, no divergence arises when f — ► 0. We have 

d T U(r, r) = f5~ l dr fi(f, r) ^ sin(2^f) exp(-2^f ) , (12) 

pr 

d 2 n(r,r) = /3- 2 9|n(f,f) ^ (4sin(27rf)exp(-47rf) - cos(2vrf) exp(-2^f)) . (13) 

Thus all terms in F[g, II] containing d T H or 9 2 II give rise to finite contributions to the radial 
integral. (The measure is dff 2 .) The same holds true for the term with (<9 r n) 2 since 

n(r, r) n(f, f) *3 1 + ?f => (a r n(r, r)f /T 2 ^ . (14) 



But 



a 2 n(r,r))^r 2 ^- (15) 

Thus 

f'°° d 2 Ii f ( \ f°° dr\ 

-J drr 2 sm(2g)-^ p [finite + 2-np 2 (lim sin(2<?) J J —j, (16) 

where the ~ sign indicates that the right-hand side approaches the left-hand side more than 
exponentially fast for increasing R > 1, see Proposition [2j Obviously, the integral in Eq. (|16p 
diverges logarithmically. ■ 

Remark 4. In summary, the quantity to be evaluated is 

.64vr 5 f ,_ 2 vr 2 / 5 4 + p 2 (2 + cos(2vrf)) f x a f°° df . 

i / dpp 1 — — i '-4, / dVL — / — sm(2o). (17) 

3 J (2n 2 p 2 + 1 - cos(2vrf)) 2 J r J R r 

Remark 5. gauge-invariant way The angular integration in the expression (|17[) . 

dQ— = / d(cos6) / dip — , (0 < a c < 2vr) , (18) 



'otC 

is regularized by introducing a defect /surplus angle ?/ < 1 for the azimuthal integration in ip: 
ac — ► aci?/' (lower integration limit) and ac — ► acTl' (upper integration limit). This singles 
out a unit vector he = (cos aci sinac, 0). Obviously, a rotation of he is induced by a rotation 
of the cartesian coordinates in which the transition to polar coordinates is performed. But for 
(j) £ K, this amounts to nothing but a global gauge rotation. Thus no breaking of rotational 
symmetry is introduced by the angular regularization. 

Proposition 6. Without restriction of generality the contribution from the anticaloron is also 
regularized in the x\X2-plane with angle aA (global gauge choice). Then we arrive at 

fC = Z c {5 al cos a c + S a2 sin ac)A(2iv(t + fc)) + £U {5 al cos a A + 5 a2 sin a a) A (2vr(f + f A ) ) , 

(19) 

where Hc,EU £ R, (undetermined: [angular integr. /xoo [radial integr. subject to dimensional 
smearing, see \44V ) an d — ^C^A < 1 (undetermined: modulus of caloron/ anticaloron which 
cannot be averaged over, see argument (Hi) above Prop. 2.). The function _4(27rf) in Eq. U9\) 
is given as 

ir 2 p 2 + cos(2vrf ) + 2 
(2ir 2 p 2 _ C os(2vrf ) + l) 2 ' 

The integral over p in Eq. A20\) diverges cubically for £ — > oo. 



A(2nr) = ¥f£dpp* 



lim sin(2<?(f , r, p)) 



(20) 
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Figure 1. The function A(^ L ) plotted over two periods with different values of £. For comparison the 



function 272£ 3 sin(2|p) is plotted as a dashed line. 



Proof. Routine computation using the fact that g saturates for p — > oo (see Prop. [2]). 

Theorem 1. The function A(2irf) rapidly approaches const^ x £ 3 sin(27rf) where const 
272.018. 



Proof. Since g saturates one may evaluate the integral numerically thus proving the claim. ■ 

Remark 6. Already for £ = 3 one has C ° nS c t onst° nSt °° = 0-025,and the functional dependence on 
f practically is a sine, see Fig.[TJ Since there is such a fast saturation towards a sine function the 
prefactor 272.018 x £ 3 , which is computed numerically, can be absorped into the undetermined, 
real number z>c,A in Eq. (fl9|) . In this sense the result for /C is independent of the cutoff £ for £ 
sufficiently large, see again Fig.CQ 



Theorem 2. The set K, coincides with the kernel of the linear differential operator T> = + 
(jf^J acting on an adjoint scalar field (f) with two polarizations. Thus T> is uniquely determined 
byJC. 

Proof. There are two independent 'polarizations' contained in fC which are given by the unit 
vectors he and h^. For each polarization there is an undetermined phase shift re A and an 
undetermined amplitude |Sc7,vl| 5 and each polarization is annihilated by T>. Modulo global 
gauge rotations there are thus two real parameters for each polarization of <p which span the 
solution space of T>(j) = 0, and T> is determined uniquely. ■ 
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T> is linear. Under the ansatz that |0| is spacetime independent Thm. [2] implies that the field <f> 
possesses a canonic kinetic term tr (d T (j)) 2 in its effective, euclidean Langrangian density. 

Theorem 3. The adjoint scalar field (f> is subject to the euclidean Lagrangian density 

£ = tr((d T ^) 2 + V(c/ ) 2 )) (21) 
with V{4> 2 ) = A 6 0~ 2 and A an arbitrary mass scale. Here (p^ 1 = -Ht^. 

Proof. Due to the BPS saturation of coarse-grained calorons/anticalorons no explicit temper- 
ature dependence may appear in 0's effective action. Together with the statement preceding 
Thm. [3] this implies an effective action of the form (|2ip with a yet unknown potential V(cp 2 ). 
The according Euler-Lagrange equations^ are [38] 



8 2 6 a = dV }}^} ^ 4> a (in components) ^ d 2 <p = dV ^t - <t> (in matrix form) . (22) 
d\4>\ z d(j) z 

Since 0's motion is within a plane in the three-dimensional vector space of the SU(2) Lie algebra, 
since \4>\ is independent of spacetime, and since 0's phase (ft is of period unity, see Thm.[2l one 
may, without restriction of generality (global gauge choice), write the solution to Eq. (|22j) as 

^ = 2|0|ti exp(±— t Z T). (23) 

BPS saturation, or equivalently, the vanishing of the euclidean energy density and Eq. (|23|) imply 



W \j) ~V(m = 0. (24) 
On the other hand, comparing d 2 (f) + (^jf^ = 0, see Thm.[2j with Eq. ([22]) . we have 

Together, Eqs. ((Ml) and fl25]) yield 

^(I0| 2 ) _ni0P) (9 ,, 
d\<P\ 2 \cf>\ 2 ■ { ° j 

Eq. (|26p is a first-order differential equation whose solution is 

A 6 

W! 2 ) = W , (27) 

where A denotes an arbitrary mass scale (the Yang-Mills scale). ■ 

Corollary 1. The modulus of the field (j) is given as = \J and hence, modulo a global 
change of gauge, 



A 3 B . 4m , . 

-^-h exp(±— t 3 r). (28) 



6 Compare with discussion in the Introduction: The search for configurations <f> with constant modulus and, 
up to global gauge transformations, pre-determined phase (j> leads to the existence of a unique and consistent (no 
contradiction to saturation of </>) value of \<f>\ and thus to the according maximal resolution. 

7 Notice our notational convention: V is either a scalar-valued function of its scalar-valued argument or a 
matrix- valued function of its matrix- valued argument. In both cases the functional dependence is identical. 
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Proof. Substitute Eq. fl27J into Eq. (|2l|) . solve for \4>\ and substitute the result into Eq. (J23]). ■ 

Remark 7. The field <fi represents coarse-grained nonpropagating, noninteracting, BPS satu- 
rated field configurations of topological charge modulus \Q\ = 1 and trivial holonomy. Thus it 
should itself not propagate. This is explicitly checked by computing the mass Mgd, of potential 
fluctuations Sep about the configuration in Eq. ([28]) as 



r2 - &v 



ML = 2 



= 48vr 2 T 2 = 12A 3 U| 2 , (29) 

2tt 



where A = 2 ^-. Since A is considerably larger than unity in the deconfining phase, see below 
where it is derived that A > A c = 13.87, and since the scale \cj)\ represents the maximal resolution 
(off-shellness) of any fluctuation after coarse-graining^] we conclude that the field 4> does not 
fluctuate: neither thermally nor quantum mechanically. The field thus represents a spatially 
homogeneous background for the dynamics of the coarse-grained, propagating gauge field (sector 
with Q = 0). 

Notice that with A c = 13.87 one obtains > 8.221 x (j^Y * , (A > A c ). But for f = 

/ \ 3/2 

8.221 x ( — ) the exponentially suppressed term below Eq. ([9"]) is a correction of less than 

/ \ 3/2 

one in 10 22 ! At the same time, setting £ = 8.221 x (— J in Eq. (|20p . one is deep inside 

the saturation regime for the set /C, see Fig.[TJ Thus, with a maximal resolution in the 
effective theory (corresponding to a length scale \cp\~ 1 up to which short-distance fluctuations in 
the fundamental fields are coarse-grained over to derive the effective theory) the infinite- volume 
limit used to derive K, and in turn the differential operator T> is extremely well approximated. 

According to Rem. [7] the configuration in Eq. (|28p is not altered by interactions with the 
gauge fields in the sector with Q = 0. Thus the nonperturbative emergence of the scale A is 
not influenced by this sector. Compare this with the situation in perturbation theory at T = 
where A is the pole position for the evolution of the fundamental gauge coupling g in the sector 
with Q = 0. There, the value of A is dictated by the value of g at a given resolution, and 
two assumptions enter. First, one assumes properties of the perturbative expansion that are 
sufficiently close to those of an asymptotic series to justify the low-order truncation of the beta 
function. Second, one assumes that close to the pole the perturbative prediction for the evolution 
of g can smoothly be extrapolated to the regime where g> 1. For an extended discussion see 



For the case of SU(3) the field eft winds in each of the three (nearly independent) SU(2) 
subalgebras for a third of the period /3, for details see [7J. 

After spatial coarse-graining the effective Lagrangian density, subject to a maximal resolution 
\4>\ for propagating gauge fields, is given as 

£ cff [a M ] = tr Qg^G'- + (D^) 2 + , (30) 

where = d^dy — d u a^ — ie[a^,a u ] = G a ^ u t a , a M = a^t a is the coarse-grained, propagating 
gauge field in the sector with Q = 0, D^cf) = 3^0 — ie[a M , 0], and e is the effective gauge coupling. 

Why is this statment true? The form of the term ^ G 2 is as in the fundamental theorjjf] due 



8 Notice that only the case of pure thermodynamics is discussed here. If an external probe is applied to the 
thermal system then an additional momentum scale enters, and \ (f>\ no longer represents the only scale of resolution. 

9 Notice that in contrast to Eq. |TJ we have not absorped the coupling into the gauge field. A gauge transfor- 
mation acting on <f> and a M thus reads: 4> —* Q<pD.' and a M — > na M n f + -fid^t , Q S fund(SU(2)) or fund(SU(3)). 
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to perturbative renormalizability [TQl CEH 121 [13] , and the only gauge- invariant way to couple 
the Lagrangian density of Eq. (|21|) to the coarse-grained sector with Q = 0, which itself cannot 
generate any composite of the field a M , is to do the replacement d T <p — ► D^ip. When we say 
'the only gauge-invariant way' we only consider local effective vertices of eft with which do 
not involve three or more external legs of the latter. (Nonlocal interactions always can be 
expanded into a local series involving powers of covariant derivatives.) Such vertices would 
mediate energy-momentum exchange from the sector Q=0 to \Q\ = 1 after coarse-graining. 
This, however, would contradict the BPS nature of the field <j> and thus its very existence. But 
(f)'s existence has just been established. Therefore no vertices involving the field <p an d other 
than the mass operator contained in tr (D^(j)) 2 (unitary gauge, no energy-momentum transfer 




but massiveness of off-Cartan modes after infinite resummation of mass insertion) exist. 

The fluctuating field is integrated out loop expanding the logarithm of the partition func- 
tion about the free quasiparticle situation. This loop expansion is nontrivial due to the term 
ie\a^,a v ] in G^ u which leads to the occurence of three- vertices and four- vertices. The momen- 
tum transfer in these vertices is subject to constraints imposed by the existence of the maximal 
resolution \<p\. The evolution of the effective coupling e is determined by the invariance of Leg- 
endre transformations between thermodynamic quantities under the applied coarse-graining up 
to a given loop order, see below. 

Apart from (small) radiative corrections and modulo global gauge transformations the full 
ground state of the effective theory in the deconhning phase (taking into account the inter- 
actions between and fundamental radiative modiheations of calorons and anticalorons) is given 
by the configuration in Eq. if28|) and the pure-gauge configuration a gs = T<Vl|§ 

The following Euler-Lagrange equation for is implied by C eS in Eq. f)30|) : 



But Eq. is solved by <j) and a gs by virtue of G^ v [a^] = D u [a^}4> = 0. 

Due to G^ u [af] = D^[a| s ]</> = the ground-state associated Lagrangian density is given as 



within calorons and anticalorons lift the energy density p gs of the ground state from zero in the 
case of BPS saturation to p BB = 47rA 3 T. The fact that the ground-state pressure P SB is negative, 
P ss = — p gs , is microscopically explained by small-holonomy calorons and anticalorons having 
their BPS magnetic monopoles-antimonopole constituents [171 [181 [Ml US [261 EZ] attract one 
another under the influence of radiative corrections, for a detailed discussion which is based on 
the work in [32], see 0H2]- The case of the 'excitation' of a large holonomy, which according to 
[42] leads to the dissociation of the associated caloron/anticaloron and hence to the liberation 
of a screened magnetic monopole and its antimonopole, is extremely rare [7J. Deviations from 
the equation of state P gs = — p es , which are due to those nonrelativistic, screened magnetic 
monopoles and antimonopoles, are described in part by the radiative corrections to the total 
pressure and energy density in the effective theory. 

Proposition 7. In the effective theory the winding gauge, where <f> is as in Eq. \28\) and a 9 ^ is 
as above, and the unitary gauge, where (ft = 2 \cf>\ ts, a 9 * = 0, are connected by a singular but 
admissible periodic gauge transformation. Under this gauge transformation the Polyakov loop 
P[a h9 \ is transformed from P = — 1 to P = 1 which points out the electric Z2 degeneracy of the 
ground state and thus deconfinement. 

Proof. Since <f> — * Cl(r)(j)C^ (t) under the gauge transformation it is easily checked that f2(r) is 





(31) 




(T = (3 1 ). Thus interactions between and radiative corrections 



given as 




(32) 
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where f2(r) = exp[±27ri^3], Z(r) = ^20(r — f ) — lj I2, and f2 gl = expfi?^]- The function 6 
is defined as 

(x < 0) , 

Q(x) = { A (x = 0), (33) 




Thus J)(t) is periodic but not smooth. The periodicity of fluctuations 5a^, however, is not 
affected by this gauge transformation. Namely, writing = + <5a^, we have 

= n„ f ns^o* + j5(t - ^)z) nj, = o gl o 5a M (n^n)* . (34) 

Now fi gl f2(r = 0) = — O g iO(r = Thus the periodicity of the fluctuation 5a^ is unaffected by 
the gauge transformation induced by fi(r) (admissibility of this change of gauge). To show the 
claimed transformation of the Polyakov loop on a^ s is trivial. ■ 

One can easily show that for SU(3) the Polyakov loop P[a^ g ] forms a three-dimensional represen- 
tation of the electric center symmetry Z3, see [TJ. Thus also for SU(3) the deconfining property 
of the thermal ground state follows. 



3.2 Thermal quasiparticle excitations 

In this section we obtain the tree-level mass spectrum for emergent thermal quasiparticles in 
the effective theory, and we derive the evolution of the effective gauge coupling e. Next, we give 
analytic expressions for the temperature dependence of thermodynamic quantities on the level of 
free quasiparticle fluctuations. We also comment on the trace anomaly of the energy-momentum 
tensor. 

We refer to an excitation, which possesses a temperature-dependent mass on tree-level in the 
effective theory, as a thermal quasiparticle. 

Proposition 8. In the effective theory dynamical gauge symmetry breaking SU(2)—>U(1) is 
manifested for the sector with Q = by quasiparticle masses m a . One has 

m 2 a = -2e 2 tr[<i>,t a }[<P,t a }. (35) 
Thus, m 2 = ml = m| = 4e 2 ^^ and m% = 0. 

Proof. Since a gs = in unitary gauge formula (|35p can be read off from the Lagrangian density 
£ cff in Eq. ([30]) . and since <f> = 2 \ <p\ t% in unitary gauge the explicit expression for the mass m 
follows. ■ 

Remark 8. Imposing unitary gauge, with gauge condition = 2 \<f)\ £3, a^ 8 = 0, and in addi- 
tion Coulomb gauge for the unbroken U(l) subgroup, with gauge condition <9ja? = 0, yields a 
completely fixed gauge if the real- valued gauge function 6 in Q 3 = exp(z6>i 3 ) vanishes at spatial 
infinity. This gauge is physical because it exhibits the quasiparticle mass spectrum, the physical 
number of polarizations - three for a = 1, 2 and two for a = 3 -, and the transversality of the 
gauge field associated with the unbroken subgroup U(l). 

Remark 9. For SU(3) the unbroken subgroup is U(l) 2 and six out of eight independent direc- 
tions in the SU(3) Lie algebra acquire mass, for details see [7JIH]- 
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Notice that the number of degrees of freedom before coarse-graining matches those after coarse 
graining. Namely, for SU(2) one has three species of propagating gauge fields (Q = sector) 
times two polarizations each plus two species of charge-one scalar magnetic monopoles (|Q| = 1 
sectoJ^l) before coarse-graining and two species of massive gauge fields times three polarization 
each plus one species of massless gauge field times two polarizations each. In both cases one 
obtains eight degrees of freedom. For SU(3) one obtains 22 degrees of freedom before and after 
coarse- graining. 

In unitary-Coulomb gauge and on the level of free quasiparticles the real-time propagators 
1^1 of the fields a]i 2 and are given as 



D^\p) = -D, v 



+ 2ir5(p 2 -m 2 )n B (\p \/T) 



p2 _ m 2 _|_ ^0 

- + 27r5(p 2 )n B (\p \/T) 
pr + 10 



(36) 



.371 



where = (g^ - ^f) , P£° = Pft = PjP = , P% = 5 ij - p^/f, u = (1, 0, 0, 0) represents 
the four-velocity of the heat bath, and ns(x) = l/(e x — 1) denotes the Bose-Einstein distribution 
function. 

Because of the existence of a maximal resolution scale \4>\ the deviation of the momentum p^ 
in Eqs. (|36l37p from its mass shell is constrained as 



\p z \ < \fi\' z , (for a=3) , \p 2 - m 2 \ < \<j>\* , (for a=l,2) . (38) 



Conditions (|38|) fix the momentum transfer in a three-vertex by momentum conservation. 
The following conditions fix the momentum tranfer in a four-vertex: 

\(pi+P2) 2 \ < \<t>\ 2 , (schannel); |(p3-Pi) 2 | < |^| 2 , (t channel) ; | {p2-Pz) 2 \ < |0| 2 , (u channel) . 

(39) 

These conditional follow from the fact that massless intermediate modes in the fundamental 
theory, which do not exist in the effective theory but dress the four-vertex such that the latter 
appears to be local, may not resolve distances smaller than |49|. 

On the one-loop level in the effective theory (gas of noninteracting thermal quasiparticles and 
massless excitations) the contribution AV of quantum fluctuations (arising from terms without 
the factor n b in Eqs. ft 3613 7]) ) is negligibly small. 

We estimate AV^ by the contribution of the massless mode appropriately weighted by the 
number of polarizations of all fields a^,a^, and oj^: 



1 ]™ dpp* log (a) = ±- 2 = A 2 



\AV\<— I dpp 3 log(^) =-^ = —V. (40) 



Since A is considerably larger than unity AV can safely be neglected. 



10 According to the definition of the set K, in Sec. 3.1 topologically nontrivial field configurations with \Q\ — 1 
only contribute to the thermal ground state. As a consequence, only magnetic monopoles of magnetic charge 
modulus unity occur as their constituents. 

11 An analytic continuation r — > —it,t and r real, is performed, see [43| . 

12 In calculating radiative corrections we have checked that under \4>\ — > £|</>| in (|38p and (|39|) . where £ is of 
order unity, the results are remarkably stable. 
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Theorem 4. On the one-loop level the evolution of the effective coupling e is determined by the 
following first-order ordinary differential equation: 



24A 4 a D(2a) 



(2vr)6 i + 



24A 3 a 2 

(27T) 6 



D(2a) 



(41) 



The evolution governed 



where a = = 2ire\ 3//2 and D(y) = f?° dx —?== , = — . 

2T yy/ J0 yfx 2 +y 2 exp(^/x 2 +y 2 )-l 

by Eq. {J^ 1\ ) possesses two fixed points: a = and a = oo. The latter is associated with a 
critical temperature A c of value A c = 13.87. An attractor to the evolution exists. It is given 
as a(A) = 4\/27r 2 A~ 3 / 2 for A » A c and a(A) oc — log(A — A c ) for A \ A c . The trace of the 
energy-momentum tensor 6^ u grows as = p 1m P = 12 A for A S> A c . 

Proof. The Legendre transformation p = T^ — P between total energy density p and total 
pressure P, which follows from the (ultraviolet regularized) partition function formulated in 
terms of fundamental fields, needs to be honored in the effective theory. Since there are implicit 
temperature dependences in the parameters \4>\ and e of the effective theory for the coarse- 
grained fluctuations 5a^ the derivatives w.r.t. temperature of these parameters ought to cancel 
one another. A necessary and sufficient condition for this to take place is d m P = 0. Because one 
may neglect the quantum part, on the one- loop level P and p are determined by the thermal 
parts of the propagators in Eqs. (|36l) and (|37l) (terms with the factor n#). They are given as 
follows: 



P(A) 



-A' 



4/ 2A 4 
1 (2vr)6 



[2P(0) + 6P(2a)] + 2A 



(42) 



p(X) = A' 



( 2A 4 
I (2vr) f 



[2p(0) + 6p(2a)] + 2A 



where P(y) = / °° dxx 2 log 



1 — exp(— \J x 2 + y 2 ) and p(y) 



f^dxx 2 



\/x 2 +y 2 



(43) 



The 



exp(^/ x 2 +y 2 ) — 1 

evolution equation (|4ip follows by applying d^ a ^P to the expression in Eq. (|42p and setting 
the result equal to zero. The right-hand side of the evolution equation ([4T]) indeed vanishes at 
a = 00 (essential zero due to the exponential in the integrand for the function D(2a)) and at 

2 

a = (algebraic zero since -D(O) exists: -D(O) = ^-). Since the right-hand side of Eq. (f4Tl) is 
negative definite this equation is equivalent to 



1 



24A 3 



X^ + a)aD(2a). 



For a <^L 1 the Taylor expansion of the function D(2a) can be truncated at zeroth ordeiri. 
simplifies Eq. (|44p as 



(44) 



This 



A 3 
(2nY 



da 

x dx +a 



(45) 



and the solution, subject to the initial condition a(Aj 



ai <C 1 is 



a 



(A) = 4^/2tt 2 A- 3 / 2 ( 1 



A 



Q»A 3 
32vr 4 



1/2 



(46) 



In [50] it was shown that albeit the coefficients in the Taylor expansion of D(y) about y = diverge for orders 
larger or equal than quadratic this formal series can be resummed to a smooth function in y. For this process the 
zeroth-order coefficient serves as a boundary condition and thus is relevant. 
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Thus for A -C Aj the function a(A) runs into the attractor a(A) = 4\/2vT 2 A _3 / 2 . Since a = 
= 27reA~ 3 / 2 there is a plateau e = \/87r in this regime. Because the attractor increases 
with decreasing A the condition a <C 1 will be violated at small temperatures. The estimate 
14.61 > A c is obtained by setting the attractor equal to unity. Since the true solution in 
this regime will continue to grow with decreasing A (negative definiteness of right-hand side 
of Eq. (I4ip ) the right-hand side of Eq. (1411) will be exponentially suppressed. This verifies the 
behavior a(A) oc — log(A — A c ) for A \ A c and implies a logarithmic singularity at A c also for 
e(A). Numerically, one obtains A c = 13.87. 

As for the (dimensionless) quantity let us consider the following function h(X) [51] : 

h(A)s _£W^W. (47) 

Expanding up to quadratic order in a, one has h(X) = 1 + \^ 2 }y! ■ Substituting the attractor 

a(A) = 4\/2vr 2 A~ 3 / 2 , one arrives at h(X) = | , (A > A c ), and the claim follows by recalling that 
pgs = -4vrA 3 T. " ■ 



The existence of an attractor to the evolution of e with temperature signals the decoupling of 
the high-temperature initial situation from the low-temperature physics. According to Eq. (|46p 
this decoupling takes place in a power-like fashion in contrast to perturbation theory. Also, 
the Yang-Mills scale A is not determined by the initial value a\ = 2irei\ i (as it would in 
perturbation theory), but it arises as a purely nonperturbative integration constant owing to 
the nontrivial thermal ground state, see Thm.[3l In a physics model the initial temperature 
Ti = 4^ is naturally given by the scale where the assumption of a smooth spacetime manifold 
supporting the Yang-Mills theory breaks down. According to present consensus this scale is 
the Planck mass. The constancy of e for a -C 1 signals that the magnetic charge g = ^ of a 
screened magnetic monopoles, liberated by a dissociating caloron/anticaloron of large holonomy, 
is conserved during most of the evolution. For A \ A c both the magnetic charge g and the mass 
M mon ~ 47r e T of a screened magnetic monopole vanish, see Rem.[TJ 

In a completely analogous way one obtains for SU(3) the evolution equation 

12A «(a) + 20(2a) 

(2») 8 l + Hg£(D(a)+2D(2a)) 

The attractor for a <C 1 reads [51] a(A) = -^7r 2 A~ 3 / 2 , and the plateau value is e = t^tt- The 
numerical value for A c is A c = 9.475, and one obtains W = 24vrA 3 T = 12 A A 4 for A > A c . 



3.3 Radiative corrections 

Here we discuss the systematic computation of radiative corrections in terms of a respective 
loop expansion within the effective theory for the deconfining phase. This loop expansion has 
little resemblence with its perturbative counterpart. Each nonvanishing diagram is infrared 
and ultraviolet finite. While the former property is due to the nonpertburbative emergence 
of (quasiparticle-)mass on tree level (adjoint Higgs mechanism) the ultraviolet finiteness follows 
from the existence of a scale \4>\ of maximal resolution: In a physical gauge quantum fluctuations 
are constrained to maximal hardness |0| 2 by the thermal ground state, and their action is small 
as compared to that of thermal (on-shell) modes. Frankly speaking, this is the reason for the 
rapid converence of loop expansions. We also exhibit two calculational examples: The on-shell 
one-loop polarization tensor of the massless mode and the dominant two-loop correction to the 
pressure. 
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After a euclidean rotation po — > ip A , (po and ^4 real) the second condition in ([38]) reads 
\p 2 + 4e 2 \(j)\ 2 \ < \(fi\ 2 . Since this is never true in the attractor regime because e > y/8ir for SU(2) 
and e = -^tt for SU(3) one concludes that massive quasiparticles do propagate thermally only. 
Alternatively, staying in Minkowskian signature, the four-momentum squared of a massive mode 
needs tuning to the mass squared of about one part in three-hundred for SU(2) and similarily 
for SU(3). Again, this implies that a massive quasiparticle propagates over long distances and 
thus thermalizes. Also massive quasiparticles cannot be created by quantum processes because 
this would invoke momentum tranfers of at least twice their mass. But this is about 35 times 
larger than the maximally allowed resolution in the effective theory. 

An irreducible bubble diagram is defined by the property that no one-particle reducible diagram 
for the polarization tensor is created by cutting in any possible way an internal line of the dia- 
gram. 

It was argued in [49J that, when computing bubble diagrams contributing to the thermal pres- 
sure, the resummation of polarization tensors, subject to an insertion of an irreducible bubble 
diagram, avoids the occurrence of pinch singularities (powers of delta functions) due to a re- 
laxation to finite-width spectral functions in the thermal parts of the so-obtained real-time 
propagators. Technically, the procedure for resumming polarization tensors is to compute them 
in real time subject to the constraints (J38]) and (I39p . perform an analytic continuation to imag- 
inary time in the external momentum, carry out the resummation, and finally continue back to 
real time. 

For a four momentum p = (po,p) circulating in a loop we refer to po and \p\ as the two indepen- 
dent radial loop momenta. We denote by K the number of independent radial loop momenta in 
a given irreducible bubble diagram. 

Independent hypersurfaces Hi , (i = 1, ■ ■ ■ , h < K) , in the K -dimensional space M, K are defined 
by the property that in a whole environment U of their intersection nf=i H the normal vectors 
normal vectors hi to Hi (computed anywhere on U n Hi) are linearly independent. 

Remark 10. If h = K then it follows that H^i Hi is a set of discrete points. 

Proposition 9. For an irreducible bubble diagram containing only V^-many four-vertices and 
no three-vertices the number K of independent constraints on the loop momenta is estimated as 
K>lV A . 

Proof. The number I of internal lines in such a diagram is I = 2 V4 |52j . Because of (|38p there 
are thus 2 V4 constraints on propagating momenta and according to ([39]) at least | V4 constraints 
on momentum transfers in vertices. (For this estimate pair the four-vertices in the diagram.) 
Together this gives K>lV A . ■ 

Proposition 10. For an irreducible bubble diagram containing only V^-many three-vertices and 
no four-vertices the number K of independent constraints on the loop momenta is given as 
K = IV 3 . 

Proof. The number / of internal lines in such a diagram is 

I = § V 3 [52]. Because of ([38]) there 
are thus | V3 constraints on propagating momenta. No additional constraints arise because 
the momentum transfer in the vertex, induced by two external legs, coincides by momentum 
conservation with the momentum of the third leg. Thus K = | V3. ■ 

Proposition 11. For an irreducible bubble diagram containing only V\-many four-vertices and 
no three-vertices the number K of independent radial loop momenta is given as K = 2 V4 + 2. 
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Proof. The Euler characteristic for a spherical polyhedron reads 2 = V — I + F where V is 
the number of vertices, I the number of edges, and F the number of faces. Since a connected 
bubble diagram is a spherical polyhedron with one face removed the identification of F = L-\- 1, 
where L denotes the number of loops (or left-over faces), yields L = I — V + 1. Combining this 
with I = 2Vi and using the fact that K is twice the number L proves the claim. ■ 

Proposition 12. For an irreducible bubble diagram containing only V^-many three-vertices and 
no four-vertices the number K of independent radial loop momenta is given as K = V3 + 2. 

Proof. Combining L = I — V + 1 with I = § V3 and using the fact that K is twice the number 
L proves the claim. ■ 

Corollary 2. From Props. WUCty 111 and{T^ one concludes that 

K 4 / 1 \ K 2 ( 2\ 

— < — I 1 + — ) , (four-vertices only) ; — = - I 1 + — I , (three-vertices only) . (49) 
K 7 \ V4J K 3 \ V3J 

Since any subdiagram (obtained by cutting more than one internal line) of an irreducible bubble 
diagram is again irreducible it follows that at a fixed number V of vertices with V = V4 + V3 > 2 
the ratio is minimal for V3 = and maximal for V4 = 0. 

If the constraints ([38]) and ([39]) were equations and not inequalities one would conclude from 
Cor. [2] that the intersection of the independent hypersurfaces Hi specified by them would be 
empty for a number L of loops greater than a finite number L max . From || < | ^1 + pM = 1 
we have V^ ax = 4 which by virtue of the proof to Prop.fTTI implies that -L max < 5. 

For completeness let us investigate the generalization of the Euler characteristic for a spherical 
polyhedron without any handles to the situation of nonplanar bubble diagrams. Notice that the 
latter can be considered spherical polyhedra with one face removed and a nonvanishing number 
of handles added (genus g > 0) : 

V-I + L + l = 2 — > V-I + L + l = 2-2g, (50) 

where again / is the number of internal lines, L the number of loops, and g represents the genus 
of the polyhedral surface (the number of handles). Notice that the right-hand side of the right- 
hand side equation is the full Euler-L' Huilliers characteristics. Reasoning as above but now 
based on the general situation of g > expressed by the right-hand side equation in Eqs. (i50j) . 
we arrive at 

I s K i+ ^ (1 " 29) )' {v=Vth 

I s K i+ I (i " 29) )' {v=V3) - <5i> 

According to Eqs. (|5ip the demand |£ < 1 for a compact support of the loop integrations is 
always satisfied for g > 1 since the number of vertices needs to be positive: V4 > and V3 > 0. 
Recall that at g = this is true only for V4 > 2 and V3 > 6, respectively. We thus conclude 
that bubble diagrams of a topology deviating from planarity are much more severely constrained 
than their planar counterparts. 

Because other thermodynamic quantities are related to the pressure by (successive) Legendre 
transformations this would imply their exact calculability as well. But since (|38p and (|39p are 
inequalities the associated hypersurfaces Hi are fattened, and the situation is less clear cut. 
However, for large V3 the ratio approaches the value |, which is considerably smaller than 
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log 10 G/T 2 




0.2 0.4 0.6 0.8 1 X 

Figure 2. |^-| as a function of X = ^ for A = 1.12 A c (black), A — 2 A c (dark grey), A = 3A C 
(grey), A = 4 A c (light grey), A = 20 A c (very light grey). The dashed curve is a plot of the function 
f(X) = 21og 10 X. There is screening to the left (G > 0) and antiscreening (G < 0) to the right of the 
cusps. The massless mode is strongly screened at AT-values for which log 10 > f(X) (^p- > X), that 
is, to the left of the dashed line. For A = A c the function G vanishes identically. 



unity, in a power like-way suggesting that D^i H = for sufficiently large but finite K (or V3 or 
L). Here Hi now refers to a fattened hypersurface. We thus arrive at the following conjecture: 

The loop expansion of the pressure in the effective theory for the deconhning phase terminates 
at a finite loop order. 

Remark 11. The argument presented in favor of the truth of this conjecture apply to both 
the SU(2) and the SU(3) case. In [53] we have observed that for the SU(2) case the three-loop 
irreducible bubble diagram with V4 = 2 and V3 = 0, containing two internal lines with massless 
and two internal lines with massive particles, vanishes identically. 

Example 3. Here we provide a one-loop example for a typical radiative correction: the polariza- 
tion tensor for the massless mode with p 2 = in the effective theory for deconfining SU(2) 
Yang- Mills thermodynamics Without restriction of generality one may assume that p points 
into the 3-direction. Then the only nontrivial entries in Yl^ v are IIn = H22 = G{po = \p\,P, T, A). 
One easily checks that only the tadpole diagram with the massive modes (a = 1,2) circulating 
in the loop contributes for one-shell external momentum (p 2 = 0). That is, for p 2 = and on the 
one-loop level there is no imaginary part in the screening functional G. By a (lengthy) routine 
computation, which takes into account the constraint:^ (I39p and recalls that massive modes 
propagate thermally only, one obtains the following result: 



G_ 



<% dp+ d£ dp 

-00 Jp m (X,£,\) Jtm(X,\) JO 



e 2 (A)A 



-3 



-4 + 



4e 2 (A) 



n B (2vrA- 3 /V/0 2 + ? 2 + 4e 2 (A) 



VP 2 + £ 2 +4e 2 (A) 



(52) 



14 On the two-loop level and at p 2 



G receives an imaginary contribution whose modulus, however, is strongly 
suppressed as compared to the modulus of the one- loop result, for the ratio of two- loop to one- loop contributions 
to the pressure see [3]. 

15 The three constraints in (|39[) collapse onto one constraint for this diagram. 
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Figure 3. The ratio of the dominant two-loop correction to the pressure and the ground-state subtracted 
one-loop result as a function of A. 



where the dimensionless quantities (, m ,£,M, Pm, and pu are given as 
2X 



M 



(53) 



pm(X,£,\) 
M 



7TN2 (4e 2 ^l) 2 
XJ A3 



2vr 4e 2 =F 1 
~X A 3 / 2 



£-4e 2 



(54) 



the dimensionless quantity X is defined as X = and e(A) follows from the solution a(A) of 

the evolution equation (f44j) by virtue of e(A) = A 3 / 2 . A plot of log 10 5% i n dependence of X 
for various values of A is presented in Fig.[2j By virtue of the constraints expressed in (|53p and 
(|54p it is straight-forward to show that the function G possesses an essential zero at X = 0. 



Example 4. In Fig. [3] the dependence on temperature of the ratio of the dominant two-loop 
diagram for the pressure (with V3 = 2 and V4 = 0) and the ground-state subtracted one-loop 
result is depicted. We refrain from quoting formulas and refer the reader to 



Remark 12. Combining the results of [7], [1], and [S3], the ratio of two-loop corrections to the 
pressure to the ground-state subtracted one-loop result is, depending on temperature, at most 
~ 1CT 2 and the ratio of three-loop corrections to the ground-state subtracted one-loop result at 
most ~ 2 x 10 -7 . 



4 Preconfining phase 

4.1 Thermal ground state: Interacting magnetic monopoles/antimonopoles 

The preconhning phase is the mediator between deconfinement at high temperatures and the 
low-temperature confining phase where (dual) gauge fields do not propagate. This phase occu- 
pies a very narrow region in the phase diagrams of SU(2) and SU(3) Yang-Mills thermodynamics. 
This and the fact that lattice simulations operate at a finite spatial volume and thus have dif- 
ficulty to exhaustively capture the important long-range correlations inherent to the monopole- 
antimonopole condensing ground state are the reasons why the preconfining phase has escaped 
its detection in numerical experiments \56\ \57\ 158] . We know that at A c one direction in the 
SU(2) algebra remains precisely massless and thus propagates. At the same time we know that 
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screened magnetic monopoles and antimonopoles become massless at A c and thus are prone to 
condensation. Even though the monopole-antimonopole condensate induces a dynamical break- 
ing of the (dual) U(l) gauge symmetry this does not lead to the immediate decoupling of the 
dual gauge mode. Why and at what point this happens and how the deconfining ground state 
relates to the monopole-antimonopole condensate (tunneling) when lowering the temperature is 
subject to insightful analysis. On average, the process of gradually generating an extra polar- 
ization of the dual gauge mode by tunneling transitions between the deconfining ground state, 
characterized by a short-lived and small caloron or anticaloron holonomy, and the monopole- 
antimonopole condensate, characterized by a stable and large caloron or anticaloron holonomy 
[42, 54, 55j appears to be relevant when addressing the emergence of intergalactic magnetic fields 
in applying SU(2) Yang-Mills thermodynamics to describe thermalized photon propagation, see 

mm ei EE]. 

In this section a derivation of the thermal ground state in the preconfining phase is given 
where magnetic monopoles and their antimonopoles are pairwise condensed. The process of 
condensation is extremely subtle microscopically: As temperature approches A c from above, 
the screening of a given preexisting magnetic monopole and its antimonopole (liberated by the 
dissociation of a large-holonomy caloron or anticaloron) is suddenly enhanced (logarithmic pole 
in e(A)). Although this rapidly suppresses their mass and magnetic charge it does not yet lead 
to the formation of a stable condensate. For A < A c the average caloron-anticaloron holonomy 
gradually increases with decreasing temperature (supercooling) giving rise to more frequent 
caloron and anticaloron dissociation processes. As a consequence, the magnetic-charge screening 
of a given monopole increasingly is due to alike monopoles and antimonopoles (and decreasingly 
due to annihilating monopoles and antimonopoles inside a small- holonomy caloron/anticaloron) 
in an ever more stable condensate. Although the derivation of a complex scalar field describing 
the monopole-antimonopole condensate only relies on the limit of total screening e — ► oo, which 
takes place at A c , the formation of a stable condensate is seen to occur at a slightly smaller 
temperature. 

Remark 13. The condensate of monopoles and antimonopoles starts to form at A c , where 
e = oo, and the Yang-Mills system 'forgets' about the existence of the mass scale A since the mass 
of monopoles and antimonopoles vanishes by total screening and since there is no interaction 
between them because the magnetic coupling g = vanishes. Also, since screened monopoles 
and antimonopoles are at rest w.r.t. the heat bath (previously created by dissociating large- 
holonomy calorons and anticalorons [32]) no kinetic energy of their motion exists. Thus, if, after 
an appropriate spatial coarse-graining, the condensate is described by a spatially homogeneous 
field <p, then this field by itself must be BPS saturated, that is, its energy-momentum tensor (or, 
after integrating out any dependence on space, its euclidean energy density) vanishes identically. 

The potential for the formation of a monopole-antimonopole condensate opens up for e — > 
oo. After an appropriate spatial coarse-graining and on the level of no interactions between 
monopoles and/or antimonopoles this condensate is, in the euclidean formulation, described by 
an inert, spatially homogeneous, and BPS saturated complex scalar field ip = \(p\ exp ±2iri^ 



where \ip\ = y and A is an arbitrary mass scale. Taking interactions between monopoles 
and antimonopoles into account, the thermal ground state is described by (p and the pure-gauge 
configuration a^' ss = ^5^4 1§ of the dual abelian gauge field . Here g is the magnetic coupling, 
and the ground-state energy density p gs and pressure P es are given as p gs = — P gs = ttA 3 T. 



complex scalar field (p: 

For the effective theory (j30|) in unitary gauge out of the three directions a^' 2 ' 3 in the Lie algebra 
only a 3 propagates for e — > oo, the other two gauge modes decouple because their quasiparticle 
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mass diverges. For A < A c we define the dual gauge field a® by the coarse-grained version of a^. 
Since is a free field for A < A c this coarse-graining is trivial. In particular no (local or nonlocal) 
composites of the field a® may propagate in the effective theory for the preconfining phase, and 



a D , -» a D 



+ £fify,fit under a gauge transformation with £1 GU(1). Because of the rotational 



symmetry of the thermal system the monopole-antimonopole condensate is described by a scalar 
field ip. Since tp by itself is BPS saturated it would be irrelevant to the thermodynamics of the 
preconfining phase if it was a U(l) gauge singlet (real scalar). But then the only option for 
coupling (p to a® in a gauge- invariant way is the transformation law <p — > Q*ip (complex scalar 
field). 
ip 's phase: 

The (dimensionless) phase 6 with ip = \ip\ exp[i6] is defined by the geometrically and thermally 
averaged magnetic flux F± ;th through a two-dimensional sphere, S2,r=oo, of vanishing curvature 
(infinite radius R) induced by a monopole-antimonopole system 16 ! at zero-momentum and e — ► 
oo. In accord with Rem. [13] this is the only possible definition of a dimensionless quantity which 
does not make any reference to a scale. Consider a system of a zero-momentum monopole 
and its zero-momentum antimonopole. In unitary gauge, where independently of position the 
adjoint Higgs field of the monopole or antimonopole configuration points into a fixed direction in 
the SU(2) algebra and Dirac strings compensate for the magnetic flux through a closed surface 
surrounding the monopole or antimonopole, we introduce unit vectors x m and x a for monopole 
and antimonopole, respectively. These vectors signal the direction of the Dirac strings (both 
pointing away from the respective center of charge). Let 5 = Z(x m ,x a ) and both monopole and 
antimonopole be placed on the same side of S2,R=oo- (It can easily be checked below that this 
is no restriction of generality.) Now a single monopole or a single antimonopole, whose Dirac 
string does not pierce S , 2 i _r= o, would induce a magnetic flux F± through S2,r=oo of F± = 
It is then easy to see [7] that, (geometrically) averaging over all directions of x m and x a at a 
given angle 5, the flux F±(5) of the monopole-antimonopole system is given as 



F ± (5) = ±^-- = ±™ , (0<5<n). 
2tt e e 

After screening the mass M m+a of the monopole-antimonopole system is given as [24] M, 



(55) 



m+a 



fa-. Thus, coupling this system to the heat bath, the thermally averaged flux F± >th (5) reads 



e/3 



F±A6)=4n / d 3 P 5^(p)n B ((3E{p))F ± (5). 



(56) 



where E(p) = J M^ +a +P 2 , and ub{x) 



i 



exp[x] — 1 



denotes the Bose function. Since 



lim I exp 



8tt 2 



, 1 8?H 1 
1+ 2— + 6 



8vr k 



+ ••■ 



(57) 



one finally has 



lim F±. th (5) = ±- = — , (0<5<vr). 

e^oo 7T ZTT 



(58) 



Since the angle 6, which, after eliminating any dependence on space by the inclusion of vanishing 
spatial momentum into the average thermal flux, see Eq. (|56p . has lost its original geometric 
meaning but still ought to parametrize a periodic situation, we may set ±~ = ±5 , (0<r< 



16 A monopole is correlated with its antimonopole in the sense that the former owes its existence to the latter 
and vice versa since their origin is the charge separation enabled by the strong deformation of a small- holonomy 
caloron/anticaloron. 
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(i). That is, periodicity of the flux in dependence of an angle in the monopole-antimonopole 
condensate is, after spatial coarse-graining, promoted to the periodicity in euclidean time of the 
associated, spatially homogeneous, BPS saturated field <p. Therefore, we have 



<P = M exp 



(59) 



Since \<p\ is spatially homogeneous with the same justification as for the field <j> in the deconfining 
phase the field <p is annihilated by the linear differential operatoi0 T> = <9 2 + (it) '■ 

£>ip = Vip* = , (60) 

where ip* denotes the complex conjugate of (p. 
ip's modulus: 

No explicit temperature dependence may appear in the eucliden action for the field <p> on the 
level of noninteracting monopoles and antimonopoles and e — > oo. According to Eq. (|60p and 
because of gauge invariance one may thus write 

S v = J Q drj d 3 x (^d T cp*d T <p + \v{\ip\ 2 )^j , (61) 

where 1^(M 2 ) is a to-be-determined gauge-invariant potential and \ip\ 2 = ip*<p. By virtue of 
Eq. ([59]) the Euler-Lagrange equation, which follows from the action (f6T|) . reads 



2 <?V(|y| 2 K.jawa /2»\ 2 dV(M 2 ) 

K ' f = ^f v UJ = -^wF- (62> 

On the other hand, the field (p is BPS saturated (vanishing of the euclidean energy density). 
Eq. ([6T]) and Eq. ([59]) thus implies that 



2tt\ 



2 



i^l 2 ( J ) -V(Wf)~0. (63) 



(64) 



Together, Eqs. ([62]) and ([63]) yield 

ay(|^| 2 ) _ ni^i 2 ) 

The solution to the first-order equation ([64|) reads 
A 6 

V(\<P\ 2 ) = J^, (65) 

where A is a mass scale which appears as a constant of integration. Substituting Eq. (|65p into 
Eq. <|S3D yields 



The quantity |(^| sets the scale of maximal resolution in the effective theory. An *S , 2 , J R=|y?|— 1 
separating a monopole in the interior from its antimonopole in the exterior (or vice versa) expe- 
riences the same magnetic flux as an S*2 _R=oo since in the condensate the monopole-antimonopole 



17 By a global U(l) gauge rotation a phase shift r — > t + To can be introduced (global rotation of the Dirac 
strings), and \ip\ so far is an undetermined normalization. This freedom spans a two-dimensional vector space 
which coincides with the kernel K, of T> and thus determines T> uniquely. 
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distance and their core-size is nil. Thus monopole and antimonopole cannot probe the finite 
curvature of S^i?^^- 1 and the infinite-surface limit is trivially saturated in the spatial coarse- 
graining. 
(p 's inertness: 

By virtue of Eqs. ([65]) and ([66]) one has 



6A 3 M 2 = 24 vr 2 T 2 , 



(67) 



where A = We will show below that A > 7.075. Thus the field ip neither fluctuates quantum 
mechanically nor thermally. 
full action and af?' gs : 

Since the field a 3 does not interact with itself the coarse-grained field a D obeys the same form of 
the action. Also, local U(l) gauge invariance dictates that d r — > 
action for the preconfining phase thus reads 



<9 M + iga®. The effective 



S 



P 



dr / drx 



1 A 6 



2M S 



dn a v ~ 9 v a^ . Making use of the inertness of the field <p, the Euler-Lagrange 



where G D V 

equations, which follow from the action (|68p . are given as 



(69) 



By virtue of V v tp = the pure-gauge configuration a 



Tjfid^ solves Eq. ([69]) . Inserting 



a^' ss and ip into ([68]) one reads off the ground-state energy density and pressure as p s l 
vrA 3 T. 



The Polyakov loop, evaluated on the dual gauge-held conhguration af!' gs , is unity independently 
of the choice of admissible gauge. 

Let us show this. The field <p remains periodic under ip — ► Wip (admissible change of gauge) 



if and only if Q = exp i I 27rn^ + a(x) I where n G Z, and a is a real function of space only. 



Hence aP. 



'P 

D + 45?<5u4 + 87 " Sun under S7, and thus the periodicity of a® is (trivially) 



assured. (Here j = 1,2,3.) In particular, 

,d, S s ^ , 2™ , djajx) 



2yr(n^ 1) 



^4 + 



dja(x) 
9 



(70) 



Thus in any admissible gauge the Polyakov loop P on a^' gs is unity: 
P[a°' gs ] = exp \ig dr af' gs ] = 1. 



The electric Z2 degeneracy of the ground state, which occured in the deconfining phase, no 
longer exists in the preconfining phase. Since the magnetic coupling g remains finite inside this 
phase this does, however, not imply complete confinement since the dual gauge field, albeit 
massive, still propagates. 

In a way completely analogous to SU(2) one derives for SU(3) the following effective action 
for the preconfining phase [7]: 



1=1 



P 



dr / crx 



1 



1 



G^G^ l + -(v^ m yv^ m + 



1 A 6 



2 VI I 



(71) 



Since SU(3)— >U(1) 2 in the deconfining phase there are now two independent species of mag- 
netic monopoles, the dual gauge fields, a®\, a^ 2 > an d the monopole-antimonopole condensate, 
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represented by inert complex scalar fields cpi, ip2- The magnetic coupling g and the scale A are 
universal, a®f s = a^ 2 ss = =f|^5^4, and the ground-state energy density and pressure are given 

as p BS = —P es = 2irA 3 T. The Polyakov loop, evaluated on the ground-state configurations a^ 1; 
a^ 2 , is unity in any admissible gauge also for SU(3). This shows that the electric Z3 degeneracy 
of the ground state of the deconfining phase no longer persists in the preconfining phase. 



4.2 Thermal quasiparticle excitations of the dual gauge field 

Proposition 13. In the effective theory for the preconfining phase the dynamical breaking of 
the residual gauge symmetry U(l) (for SU(2)) and U(l) 2 (for SU(3)) is manifested in terms of 
a quasiparticle mass m for the dual gauge field. One has m = g\tp\ = g\<fi\ = g\<f2 \ = oT where 
a = 2vr#A~ 3 / 2 . 

Proof. In unitary gauge, tp = \cp\ = ipi^ and a^' gs = a^f 8 = a^ 2 gs = 0, the relation m = g\ip\ = 
g\y>i,2\ fo r the mass of the fluctuations 5a® , 6a^, and Sa® 2 can be read off from ([68]) and (fTTj) . 
respectively (abelian Higgs mechanism), and a = 27rg\~ 3 / 2 then follows from Eq. (|66p and the 
definition A = ■ 

Remark 14. Thus, the excitations in the effective theory for the deconfining phase are free 
thermal quasiparticles. 

The contribution of quantum fluctuations to the thermodynamic pressure AV in the preconfining 
phase is negligible. 

For both SU(2) and SU(3) one obtains in close analogy to the deconfining phase the following 
estimate for the ratio ^j-: 



AV 



V 



£ 2S- < 72 > 



As we shall see, A > 7.075 (SU(2)) and A > 6.467 (SU(3)). Thus AV is a small correction to 
the (dominant! 18 !) tree-level result V . 

Proposition 14. The scales A (preconfining phase) and A (deconfining phase) are related as 

A= ( 4 + ^) ' A < (forSU(2)); A = (2 + ^-^j ' A, (for SU(SJ) . (73) 

Proof. At A c , where e = 00 and g = 0, the pressure P is continuous. Morever, no higher 
loop corrections to the one-loop result exist in the deconfining phase since in unitary gauge 
the fluctuations a^' 2 (SU(2)) and a^' 2 ' 4 ' 5,6 ' 7 (SU(3)) decouple at A c . Equating at A c = — A c the 
right-hand side of Eq. (|42p with the right-hand side of 



P(Ac) = -A 4 



-P(O) + 



(2vr) 6 K J 2 



(74) 



for the preconfining phase (negligible quantum part), yields the claim for SU(2). For SU(3) one 
needs to equate the right-hand sides of 

P(Ac) = " A4 {(2^ P( ° ) + 2Ac } (75) 



s The total pressure P is already negative at A c , see [7|- 



Yang-Mills thermodynamics 



29 



and 



P(A C ) 



-A' 



f 12% 
I (2tt)« 



P(0) + ^ 



(76) 



Theorem 5. The evolution of the magnetic coupling g with temperature is described by the 
first-order differential equation 



12X 4 



aD(a) 



(77) 



where a = 2irgX 



-3/2 



au- 



di the function D(y) is defined below Eq. j^i[ ). 



Proof. Because quantum contribution to the pressure P can be neglected in the effective theory 
for the preconfining phase one has for SU(2 



P(A) 



-A 



(2vr 6 V ; 2 



(78) 



where the function P(y) is defined below Eq. (|43p. The SU(3) pressure is just twice the SU(2) 
pressure. As in the deconfining phase, the invariance of the Legendre transformations between 
thermodynamic quantities under the applied coarse-graining implies for the effective theory that 
d( a T)P = 0, and for both SU(2) and SU(3) the same evolution equation ([771) follows. ■ 

Remark 15. Numerically, the initial condition for the evolution described by Eq. (|77p is g(X c ) = 
for A c = 8.478 (SU(2)) and A c = 7.376 (SU(3)). For decreasing A < A c the magnetic coupling 
g rises rapidly and runs into a logarithmic pole at X c r. g oc — log (A — A c '). Numerically, one 
has X c i = 7.075 (SU(2)) and X c > = 6.467. Taking the mass m of the dual gauge mode as an 
order parameter for the dynamical breaking of U(l) (SU(2)) and U(l) 2 (SU(3)) and postulating 
that m = K(T C — T) u for T ~ T c , where K and v are constants, one extracts mean-field critical 
exponents: v = ^. 

Remark 16. The energy density p divided by T 4 in the preconfining phase is given as 



rp4 



(A) 



A 4 



6A 4 

(27T)' 



A 



(forSU(2)); ^(A) 



A 4 



12A 4 



■p(a) + A 



(for SU(3)) 
(79) 



Remark 17. With a slight abuse of notation we refer to p(X) as the functional dependence of 
the energy density on temperature A in the preconfining phase and to p(X) as the functional 
dependence of the energy density on temperature A in the deconfining phase. Thus p(X) and 
p(X) are different functions of their arguments. 



At X c — ^A c 



Proposition 15 

temperature. One has A(A 



the energy density p exhibits a positive jump when decreasing the 



p(A e +0)~p(A c -0) 



| §o for SU(2) and A(A C 
3). 



8 TL 
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in for SU(3). 



Proof. Routine computation considering Eqs. 

Remark 18. The existence of the gap A signals that the monopole-antimonopole condensate 
only builds up gradually as temperature falls below A c . This is intuitively understandable 
because the condensation would require the influx of an infinite number of totally screened 
monopole-antimonopole pairs from infinity which costs energy. To facilitate the condensation of 
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additional monopole-antimonopole pairs (stable condensate) by total screening needs an increase 
of the average caloron/anticaloron holonomy from almost trivial to maximal. Then the pairs 
of liberated monopoles and antimonopoles screen one another, and no transport from infinity 
is needed. Although this process is hard to grasp microscopically, after spatial coarse-graining 
the critical temperature at which a stable condensate forms (defined by the property that the 
system is more likely to be preconfining than deconfining) can be determined exactly [3]. 

Remark 19. Notice that the number of degrees of freedom before coarse-graining matches those 
after coarse graining. Namely, for SU(2) one has one species of propagating gauge field times 
two polarizations plus one species of center-vortex loop, see Sec. El before coarse-graining and 
one species of massive, dual gauge field times three polarization. Thus, one obtains three degrees 
of freedom before and three degrees of freedom after coarse-graining. For SU(3) one obtains six 
degrees of freedom before and after coarse-graining. 



4.3 Supercooling 

A stable condensate of monopoles and antimonopoles exists for temperatures A with \ c > < A < A* 
where \ c > < A* < A c . 

At A c /, where g = oo, one has 4^ = 8ir 4 A" 3 (for SU(2)) and = 16vr 4 A" 3 (for SU(3)). 

c' c' 

On the other hand, continuing the energy density of the deconfining phase down to A c / = ^A c / 
and using Eqs. ([75]) yields 

1 cf 10 4 + 720^ 1 c/ 10 2 + 720^ 

The second summands in Eqs. (|80|) practically coincide with the above expressions for - . 

d 

Thus we conclude that p ^|'^ < p ^|'^ for both SU(2) and SU(3). But according to Prop. (p~5|) we 

have > ^ for both SU(2) and SU(3). Since both functions 4^ and nff- are continuous 

in the ranges A c ' < A < A c and A c / < A < A c , respectively, there is at least one intersection. 
Numerically, one shows that only a single intersection takes place, see Fig.^J For SU(2) one 
obtains the following values: A* = 12.15 or A* = 7.428. Now, cooling the system, a stable 
condensate (system is more likely to be found in preconfining than in deconfining state) starts 
to take place at A*, and the claim follows. 



Remark 20. The typical, maximal core-size R cor (X) of an instable center- vortex loop is given 
l 

9\<fi\ 



as i? cor (A) ~ ^ = -j-t [7J. For A c / < A < A* = 7.428, where a stable monopole-antimonopole 



condensate exists, we have g > 8.3 according to Eq. ([77]) . Thus ^ r ii^ < 0.12. That is, collapsing 

center-vortex loops are not resolved^!, and the monopole-antimonopole condensate appears to 
be spatially homogeneous. 

Remark 21. To describe the average effect of tunneling between the two trajectories 

and for A* < A < A c and A* < A < A c , respectively, one may think of the following 
'droplet' model. Let V C V tot be two volumina. The thermal probability density P(V,V tot ,X) 
for measuring a fraction of condensed magnetic monopoles and antimonopoles is given as 

exp [d(\)\ A V tot \ - 1 



The resolution is given by \(p\ 
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p(i)/T 4 "** 





11 12 13 14 15 X 



Figure 4. The quantities (dashed line) and ^J- (solid grey line) in the preconfining phase as 
functions of A. The black solid line is associated with in the deconfining phase. 



where d(X) = A(A) , and A is the temperature-dependent (positive) difference between 

trajectories and Notice that for d(X)X 3 V tot t 3> 1 the probability density P(V,V^ ot ,A) 
ceases to depend on V tot . In the model of Eq. ([51]) and for SU(2) the average polarization number 
N p of the U(l) gauge field calculates as 

/•Vtot _ / y y -V\ f Vtot - V 

^p(A)=/ dVP(V,V tot ,X)(3 — + 2^- =2+/ dVP(V,V tot ,X) — . (82) 

JO V 'tot 'tdt / JO 'tot 

Keeping V tot fixed, this yields lim^_^ N p = lim^— >o-^p = f • A similar model for the regime 
X c ' < X < A* shows that N p increases towards N p = 3 for A \ X c i. 



5 Confining phase 

In this section we turn to the confining phase which starts to set in at the temperature T c i where 
formerly instable, untwisted center-vortex loops become stable and massless [7] and the dual 
gauge field decouples because its mass diverges, see Rem. [15] At T c i the magnetic Z2 (for SU(2)) 
and Z3 (for SU(3)) symmetries start to be broken dynamically. Twisted or untwisted center- 
vortex loops, which are liberated during the subsequent decay of the monopole-antimonopole 
condensate, are interpreted as spin-1/2 fermions. 

The transition from the preconfining to the confining phase is genuinely nonthermal: Relying 
on the results of [59] for the number of connected bubble diagrams in a X<j) A theory one proves 
by Borel summation and analytic continuation that the pressure increasingly develops sign- 
indefinite imaginary parts as temperature is increased from zero towards T c i [60]. One also 
proves that at T = the pressure is precisely zero, see below [60] , 

In the preconfining phase closed lines of magnetic flud^l form by the collective dissociation of 
large-holonomy calorons/anticalorons. Inside their cores magnetic monopoles travel oppositely 
directed to their antimonopoles along the direction of the flux. The magnetic flux -F^o through 



20 Since there are no isolated magnetic charges in the preconfining phase (monopoles and antimonopoles are 
condensed) these flux lines cannot end and thus are closed. 
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a minimal spatial surface Ac, encircled by a close contour C, is given as [7j 

*±,o = | Q 9 (83) 

depending on whether a center- vortex flux pierces A once (i^) or whether it pierces Ac not 
at all or twice (0). Notice that F±^q does not depend on the velocity of the train of monopoles 
and antimonopoles travelling along the vortex line, for a discussion see [7]. 

The phase of the dual order-parameter for confinement, the 't Hooft-loop expectation (a complex 
field takes on discrete values. These are 0, in for SU(2) and 0, ±3f i for SU(3). 

The 't Hooft-loop operator <&(x, C) [62] is defined as the exponential of the magnetic flux of the 
dual gauge field through the minimal surface Mc spanned by an oriented and closed spatial 
curve C centered at the point x: 



<l(f , C) oc exp[ig (f> dzi af ] . (84) 



c 

Thus the expectation $ of $ changes phase under a singular gauge transformations of a® along 
the contour C mediated by a local magnetic center jump. Such a center jump is associated with 
an extra quantum of center flux, induced by a center-vortex loop piercing Mc in addition the 
multitude of vortices which had generated the finite value of <3? to begin with. The process of 
having an additional vortex pierce Mc proceeds in real time, and it is clear that the smooth 
dynamics of having $ change its phase is described by a complex scalar field even though for 
SU(2) it is sufficient for the equilibrium situation to assume $ to be a ^-charged real quantity. 
(Imagine the adiabatic limit, where no kinetic energy is associated with the piercing center- 
vortex loop, no change in energy is conveyed to $ by this process. If $ were real and the process 
were to be described by smooth dynamics then $ would have to change its sign in a smooth 
way (exhibiting a zero in the process) which is in contradiction to energy conservation.) 

Consider now a spatial circle of infinite radius S^ =OQ centered at x. The thermally averaged 
flux -F±,o;th of a system of a center-vortex loop and its flux-reversed partner at rest through 
A s r=oo is in the limit of vanishing core-size and mass (A — > X c i, g — > oo) given as 

lim F ±)0;th = 4tt / d 3 P 5^(p) n B ((3 c ,2E v (p, A c ,)) F ±>0 = \ x 3 / 2 (85) 

V IT ' 

where E V (Q, A) ~ ir is the typical mass of a single center- vortex loop at temperature A [7]. 
Notice the use of the Bose function n b for the system of two center- vortex loops of opposite-flux: 
Even though each vortex loop is interpreted as a spin- 1/2 fermion (two polarizations also in the 
case of selfintersections EI]) the system is of spin zero. The value zero in Eq. (|85|) is realized 
if A s r=s is pierced an even number of times by the center-vortex loops in the system, and the 

A 3/2 

values i- 2 ^— correspond to odd numbers of piercings. It is obvious that in the SU(2) case an 

A 3/2 

identification of i-^— takes place which is not true for SU(3). Properly normalized, the discrete 
values in Eq. (|85p are phase changes in for the creation of a single center-vortex loop. In SU(2) 
they are from to iir and from m to 0. For SU(3) the process to and ii-j- to create 

two distinct species of center-vortex loops. Since the spatial extent of a given center-vortex loop 
is unresolvable (g f oo, for discussion see [7]) and since in the condensate the distance between 
a center-vortex loop and its flux-reversed partner is zero, the vanishing-curvature situation is 
trivially saturated at finite curvature, S^ <oc . 

The process of decay of the monopole-antimonopole condensate(s) and the formation of the 



Yang-Mills thermodynamics 



33 



center-vortex condensate is described by real-time dynamics of the order-parameter subject 
to the potentials 



V(*) = — - A$ — - A$ , (for SU(2)) , (86) 




wLere A ~ 2 x /3 A (for SU(2)) and A ~ A (for Sl/(3jj. 

At the onset of center-vortex loop condensation thermal equilibrium is maintained at overall 
negative pressure. Periodic BPS saturated trajectories^] along euclidean time, describing the 
onset of vortex-loop condensation, exist for the potentials in Eqs. ([86]) and (|S7J), see [63]. As in 
the other two phases, this periodicity is due to the pole-term in the 'square-root' of V which 
endows the field <& with a winding number. (The 'superpotential' W has a branch cut, see |63j.) 
Furthermore, the potential V needs to satisfy the following requirements: (i) invariance under 
(local) magnetic center transformation only (no larger symmetry), (ii) dynamical realization of 
the latter (flux creation, negative tangential curvature for jump-like behavior in real time), (iii) 
the only minima of V are center-degenerate and at zero energy density (center- vortex loops in 
condensate do not interact and are massless), and (iv) as in the other two phases, a single mass 
scale A enters V. It is easy to check, see also [6l], that modulo U(l) invariant rescalings and 

~ K 2 _ A -2(n-i)^$)nj ; (k > 0,k = 1, 2, 3, • • • ,n G 25), 

which do increase the curvature of V at its minima, the potentials in Eqs. (|86p and (|87p are 
unique. Demanding at the onset of the condensation of center-vortex loops that the (negative) 
pressure be continuous in the euclidean formulation, the above relation between the scales A 
and A follows. 

Remark 22. Writing = |<3?| exp[i-|], one has 



d 2 e V(<5>) 



8 ( SU ( 2 )) 

18 (SU(3)), 1 j 



where $ min = ±A (for SU(2)) and $ min = A exp[^] , (k = 0, 1, 2) , (for SU(3)). Since |$| min is 
the scale of maximal resolution once $ has settled into one of its minima we conclude that the 
field $ does no longer fluctuate. This, in turn, implies that no tunneling to another minimum 
(flux creation) takes place once <3? has settled into $ min [7]. 

Naively, that is, without taking into account contact interactions between and internal exci- 
tations within (twisted) center-vortex loops, the thermodynamic SU(2) pressure is estimated by 
the following asymptotic-series representation^ 

\ 1=0 n>l I 



21 A canonic kinetic term in $'s effective action 

s = Jd\ (~ ( d ^yd^~\v 

is inherited from the effective action for the field ip: At the onset of vortex condensation thermal equilibrium 
prevails, and the vortex condensate coincides with the monopole-antimonopole condensate. 

22 Apologies for introducing the variable A twice in this paper, here with a different meaning than in Secs.[3]and 

m 
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Figure 5. The core of a selfintersection in a center vortex loop. 




Figure 6. Untwisted and twisted center- vortex loops up to n — 3 



where (5 = yf, M ~ A, A = e~ /3 ', L < oo, and a t £ TL. 

According to Rem. [22] no (naive) contribution to the pressure arises from the ground state: The 
fermionic gas has thermalized to a given temperature by the decay of the monopole-antimonopole 
condensate, and the field <J> is settled into one of the minima of the potential V. Because of 
its two polarization states (two directions of center flux for both SU(2) and SU(3) inherited 
from an untwisted progenitor center- vortex loop) any center- vortex loop with n selfintersections 
(n = 0, 1, 2, • • • ) is interpreted as a spin-1/2 fermion. Its mass is nM where M corresponds to the 
mass of a single intersection point (a %2 or a Z3 monopole or antimonopole is associated with the 
core of the flux eddy, see Fig. [5] marking the intersection - a plastic visualization of the concept 
of spin) which, in turn, is comparable to the scale A. Since there are two possible charges 
of the monopole singled out in the core of the intersection there are C n = 2 n many possible 
charge states of a center-vortex loop with n selfintersections. Topologically, the multiplicity 
N n of these solitons is known exactly up to n = 6 [65J, see Fig. Q. For n S> 1 the form 

Nn ~ \ ^2d={) a i n j n - 16 n was obtained in [59] by an analysis of the ground-state energy of the 
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anharmonic quantum mechanical oscillator (A</> 4 -theory in one dimension). Naturally, a; £ 2. 
Taking into account the spin degeneracy, the total multiplicity M n of a center-vortex loop with 
n selfintersections is given as 



M n = 2 x N n x C n n = X 2 x T x a* n\ 16 n . 



(91) 



Separating off the massless sector (single center-vortex loops) and negelecting any interaction 
and internal excitability, one has 



p - = i^~ 4 ( iio + ^ s Mn r dx x2 bg i 1 + 



\ n>l 



2 -/3Vn 2 +x 2 



1=0 n >l 



Here K2 denotes a modified Bessel function. In Eq. (|92p the first < sign holds strictly for the 
linear truncation of the expansion of the logarithm about unity, and the ~ sign indicates that 
terms of order (/3n) _1 have been neglected in the asymptotic expression for the Bessel function. 
This is relevant for studying the analyticity structure of the Borel resummed series. The second 
< sign holds because we have made use of the large-n expression for N n of Eq. ([91]) . Obviously, 
the expression in Eq. (f90|) represents an asymptotic series in A (zero radius of convergence) . Upon 
sending A — > —A in Eq. (|92p notice the formal similarity to the expansion of the ground-state 
energy of an anharmonic quantum mechanical oscillator [66, 67] (Ay? 4 -theory in one dimension) 
for which Borel summability was proven, see [65J and refs. therein. 

The fact that the (naive) partition function diverges because of an over-exponentially in 
energy rising density of states is known to be associated with a so-called Hagedorn transition 

EH. 



The asymptotic estimate in Eq. d90j) is Borel summable for A < 0. 

Notice that the case A < corresponds to an analytic continuation from positive-real values of 
$ = $1 + i$2 to complex values: /?2 = — > $2 = i 71 "- 



Let us now show the above claimed Borel summability. Defining P mass (A) = P as — 



and I A = 32A, the Borel transformation of P mass (A) is given as 
L L 



P_(A) = J>£A"n!n^ B P ^ (A) = £ *" n ^ • ( 93 ) 



1=0 n>l 1=0 n>l 



3 Apologies for introducing the variable A twice in this paper, here with a different meaning than in Sec.[4] 
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Thus, 5p mass (A) is a superposition of polylogarithms: 

L 

5 P mass W = 5>LL (|+0 (A). (94) 
1=0 

The functions Li_^3 + a(A) are real-analytic for A < 1. To perform the inverse Borel transforma- 
tion 

L poo 

P mass (A) = V a z £(A) = / dt e"* Bp^JX t) , (95) 

where 

/•oo 

Pi(X) = J dt e"* Li_( a +z )(At), (96) 
we notice the following integral representation of Li<,(z), valid for all s, z £ C [69J: 



Li,(z) = - / d« 7 ( , Z) " ; - , (97) 

V ^ 2 i c (l + u) s sin(vra) ' V ; 

where the path C is along the imaginary axis from — ioo to +ioo with an indentation to the left 
of the origin. Inserting Eq. (|97p into Eq. (|96|) for A = — |A| < and interchanging the order of 
integration, we have 

P t (A) = -ij du f^^l e"- " lo s(" X ) T(n + 2) . (98) 



Since, by Stirling's formulazj the gamma function T(u + 2) decays exponentially fast for -u — > 
±i oo , the integral over u in Eq. (|98p exists and defines the real-analytic El function Pi(X) , 
(A<0). 

The function P/(A) is analytic for a much larger range A G C Notice, however, that for 
arg(A) — > _±7r a branch cut is expected for Pi(X). Also, one shows that P/(0) = by using 
e (l+«)lo g (-A) = _T e «io g (-A) in Eq.JSB}- 

The approximate behavior $/(A) of Pi(X) is suggested^! as follows: 



®i Cx) = ^ r= i7 —: — > (99) 



Er=0 Q2r+l(lQg(-72r+lA)) 

Efio/32 S (log(-5 2s A)) 2 ^ 
where 72r+i,^2s £ R+, a2r+i>/?2s £ R and Si = Ri + 1. Numerically, one has for example 

log(-0.154A) 



$ (A) = 0.0570 



1 + 0.220(log(-0.494A)) 2 



a = 0.0212 log(-10.2A) + 0.00142(log(-0.109A)) 3 
U j l + 0.128(log(-1.09A)) 2 + 0.0544(log(-0.886A)) 4 ' 1 ' 



24 r(z) = \/2^z z "i e" z e H(z) where = E„>o (( z + n + V 2 ) lo S f 1 + IT^t) ~ converges for z G C_ 

and limui-Kx, -ff(z) = 0, see [70] . 

25 This follows from Eq. flUJ and the fact that Im \Li_^ 3 +i ) (2)] = for z < 0. 
26 One has: 

/" efee _a:E sin(log(-A» = a , / (fee" 1 * cos(log(-A» = lo ^ . A \.„ , {a > 0) . 

Jo a 2 + (log(-A)) 2 J a 2 + (log(-A)) 2 
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Since for <&i one has ^(A = 0) = due to a higher power of the logarithmic singularity in 
the numerator than in the denominator it is clear that |Im<J>j| grows slower than |Re<3?;| for 
sufficiently small, real-positive values of A increasing from zero. Also, Re<3?; is continuous across 
the branch cut. The growing importance of imaginary contaminations of the physical pressure 
with increasing temperature signals the growing deviation from genuine thermal equilibrium: A 
sign-indefinite imaginary part implies the existence of exponentially fast growing and decaying 
plasma modes and thus turbulences. 

Since the only difference for SU(3) is the occurrence of two types of center-vortex loops one 
obtains the result for the latter by simply multiplying the SU(2) result by two. 

6 Conclusion 

A detailed discussion and partial analysis of the thermodynamics of SU(2) and SU(3) Yang-Mills 
theory has been given. As for the case of SU(2) there appears to be a wealth of applications in 
particle physics [8] , cosmology [2j EJ HJ [5] , and plasma physics [60j [71] (dark energy by virtue 
of the axial anomaly, leptons and their interactions). To enable future contact with experiment 
in the case of SU(3) (strong interactions) the dynamics of electric-magnetically dual gauge-group 
factors (fractional quantum Hall effect) needs to be explored in their confining phases. 

Acknowledgments 

The author would like to acknowledge joyful collaboration with Francesco Giacosa on some 
aspects of the here-presented material. I am also grateful for excellent work conducted by 
my students. Many thanks go to Nucu Stamatescu for interesting discussions, and to Markus 
Schwarz and Francesco Giacosa for their helpful comments on the manuscript. I am indebted 
to my family, in particular to my wife Karin Thier, for their persistent help and understanding 
over the last four years. I also would like to thank Frans Klinkhamer for suggesting the present 
paper and for having the far-sight and fairness to care about my institutional survival. Finally, 
I would like to thank the Referee for a very thorough assessment of the manuscript and various 
helpful suggestions for improvement. 

References 

[1] Planck M., Ann. Phys. 4 (1901) 553. 

[2] Hofmann R., PoS JHW2005 021 (2006). 

[3] Giacosa F., Hofmann R., Eur. Phys. J. C 50 (2007) 635. 

[4] Schwarz M., Hofmann R., Giacosa F., Int. J. Mod. Phys. A 22 (2007) 1213. 

[5] Schwarz M., Hofmann R., Giacosa F., J 'HEP 0702 (2007) 091. 

[6] Szopa M., Hofmann R., Giacosa F., Schwarz M., arXiv: 0707.3020 [hep-ph]. 

[7] Hofmann R., Int. J. Mod. Phys. A 20 (2005) 4123; Erraturn-ibid. A21 (2006) 6515. 

[8] Hofmann R.,Mod. Phys. Lett. A 21 (2006) 999; Erratum-ibid. A21 (2006) 3049. 

[9] Garcia Perez M., Gonzales- Arroyo A., Montero A., van Baal P., JHEP 06 (1999) 001. 

[10] 't Hooft G., Veltman M. J. G., Nucl. Phys. B 44 (1972) 189. 

[11] 't Hooft G., Veltman M. J. G., Nucl. Phys. B 50 (1972) 318. 

[12] 't Hooft G., Nucl. Phys. B 33 (1971) 173. 

[13] 't Hooft G., Nucl. Phys. B 62 (1973) 444. 

[14] Shnir Ya. M.,Eur. Phys. Lett. 77 (2007) 21001. 

[15] Radu, E., Tchrakian, D. H., Phys. Lett. B 636 (2006) 201. 



38 



R. Hofmann 



Harrington B. J., Shepard H. K., Phys. Rev. D 17 (1978) 2122. 
Nahm W., Phys. Lett. B90 (1980) 413. 

Nahm W., Lect. Notes in Physics 201 (1984) 189, eds. G. Denaro. 

Hitchin N. J., Comm. Math. Phys. 89 (1983) 145. 

Corrigan E., Goddard P., Ann. Phys. (N. Y.) 154, (1984) 253. 

Actor A., Ann. Phys. (N.Y.) 148, (1983) 32. 

Chakrabarti A., Phys. Rev. D 35, (1987)696. 

Atiyah M. F., Hitchin N. J., Drinfeld V. G., Manin Yu. I., Phys. Lett. A 65 (1978) 185. 
Lee K.-M., Lu C.-H., Phys. Rev. D 58 (1998) 025011. 
Kraan T. C, van Baal P., Nucl. Phys. B 533 (1998) 627. 
Kraan T. C, van Baal P., Phys. Lett. B 428 (1998) 268. 
Kraan T. C, van Baal P., Phys. Lett. B 435 (1998) 389. 

Brower R. C, Chen D., Negele J., Orginos K., Tan C.-L, Nucl. Phys. Proc. Suppl. 73 (1999) 557. 

Gross D. J., Pisarski R. D., Yaffe L. G.,Rev. Mod. Phys. 53, (1981) 43. 

Ilgenfritz E.-M. et al., Phys. Rev. D 66 (2002) 074503. 

Ilgenfritz E.-M. et al, PoS (Lattice2005) (2005) 306. 

Bornyakov V. G. et al., Phys. Rev. D 76 (2007) 054505-1. 

Gattringer C, Phys. Rev. D 67 (2003) 034507. 

Gattringer C, Schaefer S., Nucl. Phys. B 654 (2003) 30. 

Bruckmann F, Nogradi D., van Baal P., Nucl. Phys. B 666, (2003) 197. 

Bruckmann F, Garcia Perez M., Nogradi D., van Baal P., Nucl. Phys. B (Proc. Suppl.) 129&130, (2004) 
727. 

Bruckmann F, Nogradi D., van Baal P., Nucl. Phys. B 645, (2002) 105. 
Harland D, J. Math. Phys. 48,(2007) 082905. 
't Hooft G., Nucl. Phys. B 79 (1974) 276. 
Polyakov A. M., JETP Lett. 20 (1974) 194. 

Prasad M. K., Sommerfield C. M., Phys. Rev. Lett. 35 (1975) 760. 
Diakonov D., Gromov N., Petrov V., Slizovskiy S., Phys. Rev. D 70 (2004) 036003. 
Kapusta J. I., Finite Temperature Field Theory., Cambridge University Press (1989). 
Herbst U., Hofmann R., |arXiv:hep-th/041 12141 

Herbst U., Diploma thesis, Universitat Heidelberg (2005), arXiv:hep-th/0506004 

Giacosa F., Hofmann R.. larXiv:0704.2526l [hep-th], 

Linde A. D., Phys. Lett. B 96 (1980) 289. 

Giacosa F., Hofmann R., Progr. Theor. Phys. 18 (2007) 759. 

Hofmann R., |arXiv:hep-th/0609033| 

Dolan L., Jackiw R., Phys. Rev. D. 9 (1974) 3320. 

Giacosa F., Hofmann R., Phys. Rev. D 76 (2007) 085022. 

Weinberg S., The quantum theory of fields, vol.1, Cambridge University Press (1995), 282 & 286. 

Kaviani D., Hofmann R., Mod. Phys. Lett. A 22 (2007) 2343. 

Gerhold P., Ilgenfritz E.-M., Muller-Preussker M., Nucl. Phys. B 760, (2007) 1. 

Diakonov D., Petrov V., Phys. Rev. D 76, (2007) 056001. 

Lucini B., Teper M., Wenger U., JEEP 0502, (2005) 033. 

Lucini B., Teper M., Wenger U., J HEP 0401, (2004) 061. 

Lucini B., Teper M., Wenger U., Phys. Lett. B 545, (2002) 197. 

Bender C. M., Wu T. T., Phys. Rev. Lett. 37 (1976) 117. 



Yang-Mills thermodynamics 



39 



[60] Hofmann R., Mod. Phys. Lett. A 22, (2007) 2657. 

[61] Reinhardt H., Nucl. Phys. B 628 (2002) 133. 

[62] 't Hooft G., Nucl. Phys. B 138 (1978) 1. 

[63] Hofmann R., Phys. Rev. D 62 (2000) 065012. 

[64] Hofmann R., arXv: |hep-th/0508"212"j 

[65] Kleinert H., Schulte-Frohlinde V., Critical Properties of (j) A -Theories, World Scientific Publishing Co. (2001). 

[66] Bender C. M., Wu T. T., Phys. Rev. 184 (1969) 1231. 

[67] Bender C. M., Wu T. T., Phys. Rev. Lett. 27 (1971) 461. 

[68] Hagedorn R., Nuovo Cim. Suppl. 3 (1965) 147. 

[69] Lewin L., Polylogarithms and associated functions, Elsevier North Holland Inc. (1981), 236, Eq. (7.193). 

[70] Freitag E., Busam R., Funktionentheorie, Springer (1995), 204 (Thm. 1.14). 

[71] Giacosa F., Hofmann R., Schwarz M., Mod. Phys. Lett. A 21 (2006) 2709. 



